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ABSTRACT 


This  report  outlines  the  development  of  a  computer- interactive, 
single  specimen  J-integral  fracture  toughness  test  method.  The  test  method 
was  developed  on  a  PDP  11/04  computer  and  MTS  810  servo-hydraulic  test  machine 
using  compact  tension  specimens  machined  from  a  NlMocr  gun  steel.  The  report 
discusses  some  of  the  theoretical  concepts  of  the  J- integral,  particularly 
those  which  are  relevant  to  the  test  method  which  involves  the  measurement  of 
area  under  the  load/load-line-displacement  curve  and  the  measurement  of  crack 
extension.  The  report  also  discusses  in  detail  some  of  the  errors  within  the 
computer-interactive,  single  specimen  J-integral  test  method  and  how  these  can 
be  reduced.  It  is  concluded  that  the  test  method  permits  the  evaluation  of 
the  J-integral  fracture  toughness  of  high  strength  steels  in  accordance  with 
the  existing  (astm)  standard. 


Accec’on  For 


h— ^ 

NTIS 


CR4&I 
Dl'iC  .  TAB 

U;;a:i!'Oii.-.ced 


□ 

□ 


By  . 

D^t  ib.ition/ 


'My  Codes 


Approved  for  Public  Release 


Mji!  or 
Di-St  ;  Sr^oidl 


POSTAL  ADDRESS:  Director,  Materials  Research  Laboratories 

P.O.  Box  50,  Ascot  Vale,  Victoria  3032,  Australia 


SECURITY  CLASSIFICATION  OF  TBIS  FACS 


UNCLASSIFIED 


DOCUMENT  CONTROL 

DATA  SEEET 

REPORT  NO. 

MRL-R-1021 

AR  NO. 

AR-Q 04—838 

REPORT  SECURITY  CLASSIFICATION 

Unclassified 

A  computer-interactive,  single  specimen  J-integral 
fracture  toughness  test 


AUTBOR(S)  COXFOAAZK  AUTHOR 

Materials  Research  Laboratories 


D.S.  Saunders 

I. A.  Burch 

P.O.  BOX  50, 

Ascot  Vale,  Victoria  3032 

REPORT  DATE 

September  1986 

TASK  NO. 

ARM  84/155 

SPONSOR 

Army 

FILE  NO. 

G6/4/8-3203 

REFERENCES 

35 

PACES 

r  62 

CLASSIFICATION/LIMITATION 

REVIEW  DATE 

CLASS!  FI  CATION/ RELEASE  AUTHORITY 

Superintendent,  MRL 
Metallurgy  Division 

SECONDARY  DISTRIBUTION 

Approved  for 

Public  Release 

ANNOUNCEMENT 

Announcement  of  this  report  is  unlimited 

XEYWOXDS 

J-lntegral  Fracture  Toughness  Gun  Barrel 

Single  specimen  Test  Computer  Interactive  Steel 


COSATZ  CROUPS  09020 


ABSTRACT 


This  report  outlines  the  development  of  a  computer-interactive, 
single  specimen  J-lntegral  fracture  toughness  test  method.  The  test  method 
was  developed  on  a  PDF  11/04  computer  and  MTS  810  servo-hydraulic  test  machine 
using  compact  tension  specimens  machined  from  a  NlMoCr  gun  steel.  The  report 
discusses  some  of  the  theoretical  concepts  of  the  J-integral,  particularly 
those  which  are  relevant  to  the  test  method  which  involves  the  measurement  of 
area  under  the  load/load-line-displacement  curve  and  the  measurement  of  crack 
extension.  The  report  also  discusses  in  detail  some  of  the  errors  within  the 
computer-interactive,  single  specimen  J-integral  test  method  and  how  these  can 
be  reduced.  It  is  concluded  that  the  test  method  permits  the  evaluation  of 
the  J-lntegral  fracture  toughness  of  high  strength  steels  in  accordance  with 
the  existing  ( ASTM)  standard. 


SECURXTX  CLASSIFICATION  OF  THIS  PAGE 


UNCLASSIFIED 


cokt-ents 


Page  No. 


1 .  INTRODUCTION  1 

% 

2.  THE  J-INTEGRAL  2 

3.  OBJECTIVES  OF  THE  DEVELOPMENT  OF  THE  COMPUTER-INTERACTIVE 

TEST  METHODS  6 

4.  J-INTEGRAL  TEST  SPECIMENS  7 

r' 

5.  COMPLIANCE  CALIBRATION  FOR  THE  ^-INTEGRAL  SPECIMENS  8 

6.  HARDWARE  FOR  THE  COMPUTER- INTERACTIVE  SINGLE  SPECIMEN  Jc  TEST  9 

7.  SOFTWARE  FOR  THE  COMPUTER-INTERACTIVE  SINGLE  SPECIMEN 

J-INTEGRAL  TEST  9 

7.1  The  control  program  on  the  PDP11/04-MTS  81  o  io 

7.2  The  data  analysis  programs  on  the  VAXll/780  and  PDPil/04  12 

8.  COMPUTER-INTERACTIVE,  SINGLE  SPECIMEN  J-INTEGRAL  TESTS  AND  A 

MULTIPLE  SPECIMEN  J-INTEGRAL  TEST  ON  A  Ni-Cr-Mo  GUN  STEEL  13 

9.  FURTHER  APPLICATION  OF  THE  COMPUTER-INTERACTIVE,  SINGLE 

SPECIMEN  J-INTEGRAL  TEST  METHOD  21 

10.  DISCUSSION  21 

11.  CONCLUSIONS  27 

12.  ACKNOWLEDGEMENTS  27 


13 


REFERENCES 


28 


- 

LIST  OF  SYMBOLS 

* 

;‘;»4 

a 

M 

crack  length 

. 

*0 

a 

initial  crack  length 

af 

a 

final  crack  length 

Aa 

a 

a£  -  aQ  -  crack  extension 

A 

a 

area  under  load/load-line-displacement  record 

b 

a 

uncracked  ligament  length 

B 

a 

specimen  thickness 

CTOD 

a 

crack  tip  opening  displacement 

E 

a 

Young's  modulus  S' 

E' 

a 

E/(l  -  r2) 

Err% 

a 

error  in  measurements 

‘•J* 

E(  8  ,n) 

a 

non-dimensional  function  of  t  and  n 

r 

r  V' 

xn 

a 

coefficient  of  strain  hardening  stress  intensity 

S  ^r  • 

J 

a 

J-integral  value 

\-“’y 

Jc 

a 

critical  value  of  J 

£* 

Jlc 

a 

critical  value  of  J  measured  in  Mode  I  opening 

K 

a 

stress  intensity  factor 

« 

H 

O 

a 

critical  value  of  K  measured  in  Mode  I  opening  and  by 
definition  is  used  as  a  measure  of  plane  strain  fracture 

toughness 

kq 

a 

critical  value  of  K  measured  under  conditions  which  cannot  be 
demonstrated  to  meet  all  criteria  for  a  valid  kIc  test. 

L 

a 

half-span  of  a  beam  In  three  point  bending 

;•  -; 

M 

a 

bending  moment 

■wM 

,  ■  ■  ■** 

.<V* 

a 

a 

work  hardening  exponent  in  the  Ramberg-Osgood  power  hardening  law 

'  • 

P 

a 

used  to  denote  load-control 

P 

a 

force 

r 

a 

distance  ahead  of  a  crack  tip  in  r,e  polar  co-ordinate  system 

ds 

- 

increment  in  the  contour  path 

r  _  .•; 

-  the  outward  fraction  vector  on  ds. 

-  the  displacement  vector  at  ds 

-  potential  energy  of  a  cracleed  body  under  load 

-  used  to  denote  displacement  control 

-  load-llne-dispiacement 

-  loading  work  per  unit  volume 

* 

-  constant  In  Ramberg-Osgood  power  hardening  law 

-  the  path  of  the  integral  which  encloses  the  crack  tip 

-  strain 

-  yield  strain  v 

-  strain  tensor  ^  - 

-  •  ,+•,«  total  angular  deflection  of  a  beam 

61  pi 

»  elastic  contribution  to  total  angular  deflection 

-  plastic  contribution  to  total  angular  deflection 

-  angular  displacement  ahead  of  a  crack  tip  in  r,»  polar  co¬ 

ordinate  system 

-  a  parameter  of  value  *  l  for  metals 
•  Poisson's  ratio 

-  stress 

-  yield  stress 

-  0.2%  proof  stress 

-  ultimate  tensile  stress 

O  +  <J 

0 . 2  UTS 

-  - - -  •  an  effective  yield  stress 

-  stress  tensor 

-  non-dimensional  function  of  •  and  n 


A  COMPUTER- INTERACTIVE ,  SINGLE  SPECIMEN 


■* 

J-INTEGRAL  FRACTURE  TOUGHNESS  TEST 


1 .  INTRODUCTION 

r-  ' 

In  the  fields  of  engineering'  and  materials  science  the  concept  of 
linear-elastic  facture  mechanics  Is  an  accepted  part  of  material  testing.  The 
use  of  plane  strain  fracture  toughness  as  a  material  property  Is  useful  in 
engineering  studies  because  It  allows  the  prediction  of  critical  crack  sizes 
where  the  stresses  In  a  body  or  structure  are  known.  The  measurement  of  plane 
strain  fracture  toughness,  Kjc,  is  covered  by  ASTM  Standard  E  399  C 1 J .  This 
test  is  relatively  simple  ana,  for  the  most  part,  easily  understood.  For 
plane  strain  fracture  toughness  to  be  measured  in  any  test  a  most  Important 
criterion  which  must  be  met  is  that  the  specimen  size  be  sufficiently  large  to 
ensure  that  plane  strain  conditions  prevail  across  most  of  the  specimen 
thickness.  This  means  that  there  is  a  minimum  specimen  size  for  each  material 
tested.  Thus  the  specimen  thickness  and  the  crack  length  must  both  be 
greater  than  2.5(kIc/oq  2)2,  where  Klc  Is  the  measured  plane  strain  fracture 
toughness  and  oQ  2  is  the  0.2%  proof  stress  of  the  material. 

In  many  Instances  materials  are  tested  from  parent  bodies  which  do 
not  allow  specimens  large  enough  for  the  measurement  of  Klc  to  be  taken  from 
them.  For  example,  some  modern,  large  calibre  gun  barrels  are  of 
sufficiently  high  toughness  that  valid  plane  strain  fracture  toughness 
specimens  cannot  be  taken  from  the  walls  of  the  gun  tube.  Tests  using  sub¬ 
sized  specimens  exhibit  unacceptably  high  plasticity  In  their  test  records  and 
hence  alternative  test  techniques  to  the  Klc  test  method  must  be  used,  one 
method  which  admits  a  degree  of  plasticity  is  the  J-integral  test  method  which 
Is  covered  by  ASTM  standard  E  813  [21.  The  test  method  predicts  a  critical 
value  of  J  at  which  the  crack  first  begins  to  extend.  By  convention  this 
critical  value  of  J  is  denoted  jj  because  the  test  is  predominantly  in  mode  I 
opening. 


Under  most  circumstances  the  measurement  of  Jlc  for  a  high  toughness 
material  (such  as  gun  barrel  steel)  is  an  involved,  multiple  specimen  test 
method  because  the  critical  value  of  J  has  to  be  inferred  from 
the  J-Aa  behaviour  obtained  from  a  number  of  specimens.  Aa  is  the  increase 


in  crack  length  associated  with  each  J  measurement,  and  will  be  discussed  in 
the  following  section.  However,  now  that  computer-interactive  test  machines 
are  readily  available  to  researchers  it  is  possible  to  make  use  of  the  single 
specimen  J-integral  test  method  which  utilizes  unloading  compliance 
measurements  from  which  the  extent  of  crack  growth  is  inferred. 


A  short  paper  covering  the  initial  development  of  the  computer- 
interactive  single  specimen  J-integral  test  method  at  MRL  has  been  published 
(31.  The  present  work  covers  the  same  general  aspects  of  the  test  method  as 
(31  but  expands  considerably  on  details  of  computer  programs  and  error 
analyses.  This  is  facilitated  by  practical  examples  of  the  J-integral  test 
using  a  gun  barrel  steel. 


THE  J-INTEGRAL 


The  development  of  the  concept  of  the  J-integral  as  a  parameter 
which  can  describe  the  fracture  toughness  of  a  material  is  based  on  the 
demonstration  that  a  path-independent  integral ,  J  (41,  can  uniquely  describe 
the  crack  tip  stress  singularity  in  a  material  which  obeys  the  Ramberg-Osgood 
power  hardening  law  of  deformation  plasticity.  This  power  hardening  law  for 
the  case  of  simple  tension  can  be  stated, 


<?->n 

t  y  oy  oy 


where 


stress 


-  strain 


-  the  yield  stress 


-  the  yield  strain,  7y  -  oy/E 


-  Young's  modulus 


a  material  constant 


-  work  hardening  exponent 


From  the  calculus  of  variations,  Rice  (5)  showed  that  J,  the  path- 
independent  integral,  could  be  formulated. 


J  -  Jj/wdy  -  T.(3u/ax)ds) 


ox* 


where  T.(9u/9x) 


the  rate  of  work  Input  from  the  stress  field  into  the  area 
enclosed  by  r. 


r  -  the  path  of  the  integral  which  encloses  the  crack  tip. 

W  -  loading  work  per  unit  volume, 

ds  -  increment  in  the  contour  path. 

T  -  the  outward  traction  vector  on  ds . 

u  -  the  displacement  vector  a£  ds. 


J  can  also  be  interpreted  as  the  change  in  the  potential  energy  of  a 
cracked  body  with  a  small  increment  in  crack  length  C61. 


Using  the  path  Independent  characteristics  of  the  J-integral,  Rice 
and  Rosengren  (71  and  Hutchinson  [8a,bJ  showed  that  the  following  equations 
for  the  stresses  and  strains  at  the  crack  tip  apply; 


o  -  ;  E - - — - — ll/n+1  E  (,  n) 

11  y  «5y  7y  V  15 


-  .  j  .n/n+i  r  , .  _  4 

«li  -  «<y  f - 1  (*,n) 

J  aoy  7y  lflr  J 


where  oy  »  mean  material  yield  stress 
7y  -  mean  material  yield  strain 

«  -  constant  in  the  Ramberg-Osgood  power  hardening  law 

n  -  work  hardening  rate  from  the  Ramberg-Osgood  power  hardening 
law. 

In  -  coefficient  of  strain  hardening  stress  intensity  which  has 
been  evaluated  numerically,  [8a, b). 

r,e  -  polar  co-ordinates. 

and  J  can  be  evaluated  for  a  boundary  removed  from  the  crack  tip. 

The  similarities  between  the  above  formulations  in  terms  of  J  and 
those  in  terms  of  K,  the  stress  intensity  factor  t9],  are  readily  apparent; 
see  also  Liu  tiol,  and  Figure  l  for  a  summary  of  the  HRR  model  and  the  linear- 
elastic  model  of  the  stress-strain  fields  at  the  crack  tip. 


we  see,  above,  that  J  can  be  used  to  describe  the  stress  and  strain 
singularities  at  the  crack  tip  and,  like  K,  J  can  be  measured  by  a  mechanical 
test,  thus  J  can  be  evaluated  physically  as  the  change  in  potential  energy,  u, 
as  the  crack  extends;  see  (ill  and  [12]  for  the  following  proof. 


J 


(|2> 

3a  p 


and  J 


(|2) 

3a  v 


(5) 


where  p  and  v  denote  load  and  displacement  control  respectively. 
For  the  linear-elastic  case  we  write: 


_  1  _  3v  .  _ 

J  -  —  P  —  and  J 
2  3a 


-± 

2  3a 


(«) 


and  for  the  non-linear  elastic  case: 


j  -  fp  dp  and  J 

Jo  3a  P 


-J 


V  ,3P. 

(— )  dv 
0  3a  v 


(7) 


If  we  deal  with  deeply  cracked  beams  we  can  write  the  above  in  terms 
of  the  deflections  of  the  beam,  where: 


®T  " 

9  +8 

el  pi 

where  8_ 

T 

-  total  deflection 

®el 

■  elastic  contribution 

to 

total 

deflection 

V 

-  plastic  contribution 

to 

total 

deflection 

When  the  ligament  is  fully  plastic,  eT  -  •  and  .‘.  «T  -  v/L  (9) 

where  v  -  back  face  displacement  of  the  beam 
L  -  1/2  the  loading  span  of  the  beam 

M  % 

thus  9  -  f( — - ,  —,  n)  (10) 

p  Bb  o 

y 


4 


thus  P 


(12) 


( 13) 


J  -  ~  fV  Pdv  (14) 

D  Jo 


we  note  that  JV  Pdv  Is  the  area.  A,  under  the  load/load-point-displacement 
curve  for  the  °beam. 

Equation  14  Is  the  fundamental  relationship  which  is  used  for  the 
evaluation  of  J  for  small  Increments  of  crack  length.  This  is  used  in  the 
present  ASTM  standard  for  the  measurement  of  Jlc  [21‘ 

There  are  two  methods  by  which  J  can  be  measured,  one  using  a 
multiple  specimen  technique  where  a  series  of  specimens  is  tested  covering  a 
range  of  crack  increments  and  the  other  using  a  single  specimen  where  the 
length  of  the  crack  is  continuously  monitored  as  it  extends  during  a  test. 

The  J-lntegral  value  is  determined  from  the  area  under  the  load/load-llne- 
displacement  record  using  equation  14  above.  This  method  is  summarized 
dlagrammatlcally  in  Figure  2.  Both  methods  lead  to  the  measurement  of 
a  J-Aa  curve.  This  curve,  in  Itself,  does  not  determine  the  critical  value  of 
J,  Jc,  where  the  crack  begins  to  extend.  This  is  done  by  linearly  back- 
extrapolating  the  J-Aa  data  to  the  "blunting  line".  The  blunting  line 
describes  the  apparent  crack  growth  which  occurs  as  the  crack  tip  blunts 
during  the  initial  development  of  the  plastic  zone.  The  blunting  line  is 
derived  from  the  relation  given  by  Rice  and  Johnson  (13)  where: 


Bb2  -  ,v  . 

- —  o  h  (7,  r#  n) 
L  y  L  E 


,  3P 
we  know  —  - 


3P  2bB  -  ,v  _J_  . 

ib  '  V  (I-  e  '  nl 


M 

and  substituting  L  -  —  in  ( l 3 )  and  (7) 


eff 


J  »  X  o 


CTOD 


(15) 


where 


l  for  metals 


o  -  effective  yield  strength  of  the  material 
CTOD  -  crack  tip  opening  displacement. 


The  effective  increment  in  the  position  of  the  crack  tip,  Aa,  can  be 
considered  to  be  about  0.5  x  CTOD  and  hence  a  simple  relation  for  this 
blunting  line  can  be  written 


J  -  2  x  a  __  x  Aa 
eff 

°0  2  +  °UTS 

where  o  „  is  often  defined  as  { — 1 — - ) 

eff  2 


The  point  where  the  linear  regression  line  through  all 
admissible  J-Aa  data  pairs  1 2 ]  intersects  the  blunting  line  is  considered  to 
be  the  critical  value  of  J  at  which  the  crack  first  extends. 

What  constitutes  admissible  J-Aa  pairs  is  discussed  in  detail  in 


The  computer-interactive,  single  specimen  J-integral  test  method 
described  in  the  following  sections  permits  the  measurement  of  Jlc  in 
accordance  with  the  ASTM  standard  1 2  3 . 


3.  THE  OBJECTIVES  OF  THE  DEVELOPMENT  OF  THE  COMPUTER- INTERACT IVE , 


SINGLE  SPECIMEN  J-INTEGRAL  TEST  METHOD 


When  all  refinements  are  completed,  (and  this  can  only  be  achieved 
by  the  continued  use  of  the  test  method  developed  here),  the  test  method 
should  permit  the  determination  of  Jlc,  l.e.  fracture  toughness,  on  a  routine 
basis  and  in  accordance  with  ASTM  E  813  (2). 

The  objectives  of  the  development  of  the  computer-interactive  single 
specimen  J-integral  test  method  are: 


*  The  ASTM  Standard,  (21,  suggests  that  "when  estimating  ofiff  influences 
of  testing  conditions  such  as  loading  rate  and  temperature  should  be 
considered" . 


( 1 )  Accurately  determine  JIC' 

(2)  Document  all  relevant  calculations. 


(3)  Store  all  test  data  in  a  logical  and  retrievable  form, 

(4)  Be  flexible  in  a  mathematical  sense  to  allow  compliance 
calibrations,  Young's  modulus,  Merkle-Corten  coefficient  and  other 
test  parameters  to  be  inserted  into  the  program  or  readily  changed, 

(5)  Ability  to  re-examine  test  data. 


Due  to  the  limited  transient  programming  area  ( TPA)  of  32  K  words 
available  in  the  PDPll/04  computer  at  the  time  of  the  development  of  the 
program  for  the  control  of  the  J-integral  test,  the  program  may  be  regarded  as 
incomplete  because  it  does  not  undertake  the  linear  regression  through 
the  J-A  a  test  data  to  the  blunting  line  for  the  calculation  of  JIc.  This  can 
be  done  using  a  plotting  file  created  for  this  purpose.* 


The  objectives  (l)-(4)  are  met  by  the  ct^rrent  J-test  program  for  the 
PDPll/04;  JCTS.BAS,  for  the  compact  pension  specimens,  objective  5  is  met  by 
programs  written  for  both  the  VAX1 1/780  and  the  PDPll/04  which  allow  files 
written  by  the  PDPll/04  under  RTll  version  4.1  to  be  treated  as  data  files  for 
what  is  effectively  another  J-integral  test.  This  program,  JRERUN.BAS,  re¬ 
examines  the  existing  data  but  allows  the  insertion  of  accurate  specimen 
dimensions  and  crack  lengths  and  the  use  of  alternative  compliance 
calibrations  in  the  calculations.  Estimates  of  the  accuracy  of  the  linear 
regression  analysis  to  infer  crack-length  are  made  by  this  program. 

The  programs  are  listed  in  Appendix  1  and  details  of  these  programs 
are  given  later  in  this  report. 


4.  THE  J-INTEGRAL  TEST  SPECIMENS 


The  present  work  was  mainly  undertaken  for  the  "compact  tension" 
specimen  design  and  although  the  test  method  is  readily  applicable  to  beam 
specimens,  it  is  not  dealt  with  in  detail  in  this  report. 


This  feature  will  be  incorporated  in  the  program  when  the  computer  Is 
upgraded  to  a  PDPll/34  which  ha3  a  larger  TPA. 


Several  specimen  designs  have  been  proposed  for  the  compact  tension 
specimen  for  J-integral  testing  (141.  fhe  primary  requirement  from  the 
specimen  is  the  need  for  the  measurement  of  the  load-line-displacement,  see 
equation  14  above. 


For  J-integral  testing  the  present  worK  used  a  modified  compact 
tension  specimen  which  conforms  to  the  ASTM  E  399  specimen  geometry  except 
that  it  has  an  Integral  icnife  edge  machined  at  the  load-line,  see  Figure  3. 
This  specimen  is  similar  to  that  proposed  by  Clarice  et  al  C141,  except  that 
the  pin  holes  are  closer  together.  The  advantage  of  utilizing  specimens  which 
are  essentially  common  to  the  E  399  and  E  813  fest  methods  is  that,  where  Klc 
cannot  be  measured,  the  specimens  are  readily  converted  to  J-integral 
specimens.  The  MRL  specimen  design  would  probably  preclude  the  use  of 
clevises  with  roller  bearings,  however,  there  is  no  conclusive  evidence  that 
this  is  a  necessary  feature  of  J-integral  tests  using  the  compact  tension 
specimen. 


The  J-integral  specimen  manufactured  by  mrl' is  designated  mrl  J03, 

[151.* 
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5.  THE  COMPLIANCE  CALIBRATION1  FOR  THE  COMPACT  TENSION  J- INTEGRAL  SPECIMENS 


At  the  present  time  there  are  two  elastic  unloading  compliance 
calibrations  which  are  Immediately  applicable  to  the  mrl  J03  specimen.  The 
first  is  an  experimental  unloading  calibration  using  13  mm  thick  mrl  J03 
specimens  1161,  and  the  second  a  theoretical  compliance  calibration  developed 
from  a  finite  element  program,  PAFEC,  run  on  a  VAX  11/780  computer  Il7a,bl. 
These  compliance  calibrations  are  plotted  in  Figures  4  and  5  together  with  the 
unloading  compliance  calibrations  supplied  in  ASTM  E  813  (see  also  Saxena  and 
Hudak  t 1 8 1  and  Newman  (191).  These  calibrations  are  used  in  Section  8  of  this 
work.  Listings  of  the  compliance  data  for  compact  tension  specimen  designs 
are  given  in  Appendix  2. 


* 


This  is  a  specific  designation.  The  "03"  refers  to  an  a/w  of  0.3.  Other 
starter  crack  lengths  can  be  used. 


6.  THE  HARDWARE  FOR  THE  COMPUTER- INTERACTIVE 


SINGLE  SPECIMEN  Jlc  TEST 


The  single  specimen  J-integral  tests  were  undertaken  on  a  250  kN 
servo hydraulic  testing  machine  (MTS  810  system),  interactive  with  a  PDPii/04 
'computer .  * 


The  graph  generated  by  the  program,  J-Aa,  was  plotted  on  a  Tektronix 
4010-1  VDU  and  dumped  via  a  Tektronix  thermal  imager  while  the  analogue  output 
of  load/load-llne-displacement  was  plotted  using  a  HP  7090A  plotter.  This 
plotter  can  be  controlled  by  the  PDPll  and  hence  could  be  used  for  the  re- 
plotting  of  the  .data  as  stored  by  the  program.  More  conveniently  this  is 
handled  via  the  VAXll/780  and  appropriate  graphics  terminals. 


The  load-line-displacement  was  measured  using  an  MTS  632.02C.20  clip 
gauge.  The  locating  edges  on  the  clip  gauge  arms  were  modified  slightly  to 
accommodate  the  MRL  J03  integral  knife  edge  design.  The  linearity  and 
reproducibility  of  this  clip  gauge  were  tested  experimentally  and  the  data 
plotted  in  Figure  6.  It  was  considered  unnecessary'' to  use  any  correction 
factors  over  the  range  used  for  the  J-Antegral  test  of  an  mrl  J03  specimen. 


7.  THE  SOFTWARE  FOR  THE  COMPUTER-INTERACTIVE 
SINGLE  SPECIMEN  J-INTEGRAL  TEST 


The  software  written  for  the  single  specimen  J-integral  test  method 
and  the  data  analysis  programs  are  largely  interactive  and  self  explanatory, 
however  brief  descriptions  of  the  programs  are  given  in  Appendix  3.  some 
documentation  of  the  variables  in  the  control  program  for  the  PDP  ll/04,  where 
relevant  to  the  description  of  the  software  in  this  section,  is  given  in 
Appendix  4. 


The  programs  which  undertake  J-integral  testing  and  data  analysis 
are  in  two  separate  suites. 

(1)  The  computer-interactive  control  program  written  for  the  pdp  11/04. 

and 

(2)  The  data  analysis  programs  written  for  the  VAXll/780  and  pdpii/04 
which  utilize  the  data  generated  and  stored  on  disc  by  the  control 
program. 


This  computer  has  a  TPA  of  32K  words  but  is  to  be  upgraded  to  a  PDPll/34 
with  a  TPA  of  64K  words. 


7.1  The  control  program  on  the  PDPll/04-MTS  810 


The  language  used  for  the  control  program  is  dec  Multi-user  basic 
with  MTS  and  Tektronix  enhancements  for  the  Master  segment  Generator  and  the 
graphics  respectively  (20a,bl.  The  Master  segment  Generator  (MSG)  is  used  to 
generate,  through  a  D/A  converter,  the  offset  signal  to  drive  the  servo-valve. 

The  single  specimen  J-lntegral  test  is  conducted  in  stroke  control, 
with  the  position  of  the  ram  being  measured  by  its  displacement  transducer 
rather  than  using  the  displacement  of  the  specimen  as  measured  by  the  clip 
gauge.  This  makes  the  test  fail  safe  should  the  clip  gauge  signal  fail. 

At  the  present  time,  for  J- integral  tests  using  the  mrl  JO 3 
specimens,  ±  2047  bits  >  ±  5  mn  of  ram  travel.  This  stroke  can  be  decreased 
to  allow  an  Increase  in  the  number  of  data  points  acquired  for  any  test  using 
J03  specimens,  incrementing  the  stroke  during  the  test  results  in  a  small 
increment  in  the  load  applied  to  the  specimen,  while  displacement  in  the  load- 
line  of  the  notch  is  measured  by  the  clip  gauge  where  ±  2047  bits  -  ±  2  mm. 
For  the  MRL  J03  specimen  it  is  found  that  2  mm  movement  of  the  clip  gauge  is 
ample  for  most  J-lntegral  tests.  This,  calibration  can  be  altered  using  the 
MTS  440.21  signal  conditioner  to  increase  the  sensitivity  of  the  test. 


Figure  7. 


The  software  for  the  test  method  is  summarised  in  the  flow  diagram. 


The  single  specimen  J-integral  test  method  requires  the  collection 
of  load/load-line-displacement  data  for  each  bit-wise  increment  (or  decrement 
in  the  case  of  an  unloading)  produced  by  the  Master  segment  Generator, 

(MSG).  Each  bit-wise  increment  and  decrement  in  stroke  is  accompanied  by  the 
sequential  reading  number  and  with  each -reading  is  recorded  load,  load-line- 
displacement,  cumulative  area  under  the  load/displacement  record  and  a 
switching  Integer  (Z9),  which  has  a  value  of  1,2,3,  or  4  depending  on  the  call 
to  the  MSG,  i.e.  the  ramp  control.  These  data  are  written  sequentially  into  a 
large  virtual  file  on  disc  as  binary  numbers.  The  switch,  Z9,  is  the  integer 
used  to  control  the  loading  and  the  unloading  sequences  of  the  test  and  is 
Incorporated  in  an  "ON  GOTO"  statement  to  call  the  bit-wise  output  from  the 
MSG,  see  Figure  7.  The  ramp  control  is  shown  schematically  in  Figure  8.  The 
switch  has  an  important  secondary  function  in  the  JRERUN . BAS  programs  on  the 
VAXll/780  and  PDF  11/04  where  this  record  is  used  to  control  the  use  of  the 
load/load-llne-dlsplacement  data  in  a  manner  analogous  to  the  control  program 
on  the  PDPll/04,  viz.,  the  switch  designates  data  which  are  to  be  treated  as 
unloading  data  for  inferring  crack  length  via  the  unloading  compliance  and  the 
normalized  compliance/normalized  cracklength  calibration  such  as  those  in 
Figures  3  and  4.  JRERUN. BAS  is  discussed  in  detail  in  section  7.2. 

The  acquisition  of  the  load/load-line-displacement  data  is 
undertaken  as  a  separate  (but  common)  subroutine  (through  a  "GOSUB"  statement) 
called  by  each  of  the  ramp  routines  set  by  the  ramp  switch,  i.e.  the 
acquisition  of  data  is  always  undertaken  in  the  same  manner.  Thus,  for  each 
sequential  bit-wise  increment  or  decrement  in  stroke  from  the  MSG,  when  the 
ram  of  the  MTS  has  stopped  moving,  ten  reading  pairs  of  load  and  load-line- 
displacement  are  taken  at  1  millisecond  intervals  and  then  averaged.  These 
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average  values  are  "written"  Into  the  record  and  used  for  the  calculation  of 
the  cumulative  area  under  the  load/displacement  curve.  Once  all  calculations 
have  been  completed  the  ramp  is  re-activated,  see  Figure  8. 

Initially  the  test  behaves  in  a  linear-elastic  manner.  A  load,  P7, 
to  which  the  specimen  can  be  loaded  and  still  remain  linear-elastic  is  used  to 
establish  that  the  area  calculation  is  correct.  This  is  done  by  comparing  the 
cumulative  area  up  to  this  elastic  load  with  a  geometric  calculation  using  the 
slope  of  the  linear-elastic  loading  line  and  the  elastic  load;  viz.  the  area 
of  a  triangle.  In  the  linear-elastic  region. the  cumulative  area  should  be  to 
within  +2%  of  the  geometric  calculation. 

The  test  parameters  such  as  specimen  dimensions,  Young's  modulus, 
elastic  load  for  preliminary  calculations,  Merkle-Corten  coefficients  1 2 1 1  and 
the  load  cell  calibration  are  inserted  into  the  computer  and  also  stored  in 
the  data  file  ("name". DAT)  at  the  start  of  the  test.  A  plotting  file 
("name".PLT)  is  also  set-up  concurrently  with  the  data  file.  The  graphics  for 
the  J-Aa  plot  are  also  set-up  before  the  test.  The  test  starts  with  an 
initial  loading  in  the  linear-elastic  range,  the  limit  set  by  the  elastic  load 
(Z9-l),  and  the  area  calculation  is  checked  as  described  above.  After  a  20 
bit  increment  in  stroke  when  Z9-2  (see}  next  paragraph),  a  io%  unloading  (Z9-3) 
allows  the  Initial  crack  length  <aQ)  to  be  established  from  the  unloading 
compliance  calibration  in  the  program.*  To  produce  crack  growth  in  a 
controlled  way  the  bit-wise  increment  in  stroke  (Z9-2)  is  limited  to  a  pre-set 
number  of  bits  from  the  MSG.  This  can  be  changed  readily  within  the  program 
but  for  a  compact  tension  specimen  geometry  such  as  the  mrl  J03  it  has  been 
found  that  20  bits  total  output  from  the  MSG  for  each  loading  ramp  is 
sufficient  to  produce  measurable  crack  extension.  As  the  stroke  is  increased 
the  cumulative  area  under  the  load/displacement  curve  is  "written"  into  the 
data  and  plotting  files  (VFl  and  VF2  respectively) .  On  attaining  the  total 
number  of  bits  of  stroke  increment  the  specimen  is  unloaded  (Z9-3)  to  90%  of 
the  load  in  the  preceding  loading  ramp  and,  using  a  linear  regression  analysis 
on  the  data,  the  unloading  compliance  is  calculated  and  normalized  with 
respect  to  E  (Young's  modulus)  and  B  (specimen  thickness).  No  area  under  the 
curve  is  calculated  in  this  section  of  the  program.  Through  the  compliance 
calibration  in  the  program  (see  Section  5),  crack  length  is  inferred  and  the 
corresponding  J-lntegral  value  calculated  using  equation  14.  A  summary  of 
these  results  13  "written"  into  the  plotting  file  (VF2).  The  j-integral  test 
then  continues,  firstly  by  restoring  the  specimen  to  the  displacement  at  the 
top  of  the  unloading  ( Z9— 4 )  and  then  by  bit-wise  incrementing  the  stroke  a 
further  20  bits  (Z9-2)  to  promote  further  crack  growth.  The  determination  of 
the  area  under  the  load/displacement  curve  is  resumed  at  the  point  in  the 
program  when  Z9«2. 


It  has  been  found  that  aQ  is  often  better  estimated  if  this  first  10% 
unload  is  used  for  the  calculation  rather  than  the  initial  elastic  loading 
as  specified  in  (2).  This  unloading  must  occur  in  the  linear-elastic 
regime. 


The  data  "written"  into  the  virtual  files  (vfi  and  VF2)  can  be 
retrieved  using  the  program  ACCESS. BAS.*  All  details  of  the  test  Including 
scaling  factors  are  included  in  the  printout  of  data  from  the  PDPll/04  virtual 
files. 

A  typical  analogue  output  from  a  single  specimen  J-lntegral  test  is 
"given  in  Figure  9;  this  will  be  dealt  with  in  detail  in  section  8  when  the 
test  results  from  gun  steel  are  discussed. 

7.2  The  data  analysis  programs  on  the  VAXl 1/780  and  PDPll/04 

The  maiin  program  in  this  suite  of  programs  for  the  VAXll/780  and 
PDPll/04  is  for  the  re-analysis  of  the  test  data  for  the  purpose  of  arithmetic 
checlcs  and  the  use  of  alternative  calibrations  or  test  details  such  as 
specimen  dimensions  etc.  To  use  this  program  on  the  VAX  the  data  file  used  by 
the  JRERUN.BAS  program  must  be  in  ASCII  code.  The  data  file  created  on  the 
PDPll/04  and  loaded  into  the  VAXll/780  is  in  binary  numbers  and  must  therefore 
be  re-written  in  ASCII  code.  This  is  achieved  by  a  small  program  to  create  an 
ASCII  code  output  file  which  can  be  utilised  by  the  data  anlaysis  programs 
(see  Appendix  2  for  details).  In  the  tase  of  JRERUN.BAS  for  the  PDP  11/04  the 
virtual  files  are  read  directly.  In  all  other  respects  the  programs  are 
similar. 

The  data  analysis  program,  JRERUN . BAS  requires  specimen  details  to 
be  input  again  (should  errors  be  discovered  in  the  original  measurement  of  the 
specimen)  ahd  also,  with  the  accurate  measurement  of  initial  craclc  length  now 
that  the  test  specimen  has  been  broken  open  after  heat  tinting,  a  correct 
value  of  the  Merkle-corten  coefficient  can  be  inserted  for  compact  tension  J 
tests.  Similarly  Young's  modulus  can  be  changed  for  the  re-analysis  of  the 
test  data.  Most  importantly,  however,  alternative  compliance  calibrations  can 
be  used  if  they  are  available  for  the  specimen  geometry  tested.  This  should 
allow  the  eventual  attainment  of  better  correspondence  between  measured  crack 
length  values  and  those  Inferred  by  the  unloading  compliance  technique.  The 
standard,  ASTM  E  813  (21,  requires  "that  the  final  crack  extension  value,  as 
predicted  by  single-specimen  techniques,  must  agree  with  the  averaged  heat 
tint  values  within  15%".  Using  the  present  re-analysis  method  this 
requirement  can  usually  be  met  for  the  compact  tension  J-lntegral  specimen 
(see  section  9).  The  JRERUN.BAS  program  also  gives  the  correlation 
coefficients  for  linear  regression  analysis  of  the  unloading  data  and  also 
permits  the  examination  of  the  load/load-llne-displacement  data  pairs  to 
decide  if  any  should  be  discarded  from  the  start  of  the  unloading  to  allow  for 
non-linear  effects.  Some  computer-interactive  J-integral  test  programs  use  a 
delay  before  an  unloading  to  reduce  non-linear  effects  and  time-dependent, 
plastic  recovery  (22,23).  The  control  program  on  the  PDPll/04  overcomes  non¬ 
linear  effects  to  some  extent  by  discarding  a  fixed  number  of  data  points  for 
the  linear  regression  to  establish  the  compliance.  The  data  analysis  program 
on  the  other  hand  allows  the  number  of  points  which  are  discarded  to  be  varied 
at  each  unloading. 


The  program  prints  a  summary  of  all  calculations  together  with  all 
test  details  including  the  compliance  calibration  used  to  infer  the  crack 
extension  in  the  J-integral  test. 

Another  program,  ACCESSl.BAS,  for  the  VAX  only,  is  a  simple  utility 
to  allow  the  examination  of  any  part  of  the  load/load-llne-displacement  record 
‘(stored  as  load/load-llne-displacement  data  pairs).  This  program  "writes"  a 
simple  plotting  file  which  can  be  used  with  MRL  utility  program  GRAPH  [24]  so 
that  the  record  between  any  two  reading  numbers  can  be  replotted.  The 
ACCESS 1 .BAS/GRAPH  combination  can  be  used  to  examine  for  any  hysteresis  in  the 
unloading/ reloading  sequence  and  to  obtain  the^  slope  of  any  unloading  or 
reloading  part  of  the  record,  should  this  be  necessary  (see  Section  10). 


8.  COMPUTER- INTERACTIVE ,  SINGLE  SPECIMEN  J-INTEGRAL  TESTS  AND 
A  MULTIPLE  SPECIMEN  J-INTEGRAL  TEST  ON  A  Nl-^Cr-MO  GUN  STEEL 


The  material  used  in  this  wqrk  was  taken  as  MRL  J03  compact  tension 
specimens  in  the  transverse  orientation  from  a  4.5  inch  Naval  gun.*  This 
material  was  used  for  the  preliminary  investigation  of  the  J-integral  test 
method  (31  because  limited  fracture  toughness  data  was  available  for  this 
steel  (25,281.  The  chemical  composition  of  the  steel  is  given  in  Table  1. 

A  micrograph  of  the  steel,  with  orientation  transverse  to  the  barrel 
wall,  is  shown  in  Figure  10. 

The  mechanical  properties  of  this  gun  steel  have  been  reported  in 
earlier  work  (251  and  are  summarised  in  Table  2.  Estimates  of  fracture 
toughness,  using  ASTM  E  399  Ill,  were  made  from  the  load/displacement  records 
of  tests  on  the  steel  (281;  these  are  reproduced  in  reference  25.  It  should 
be  noted,  however,  that  the  fracture  toughness  values  reported  in  Table  2  are 
designated  as  Kg  values  because  the  fracture  toughness  tests  were  invalid 
(261.  For  the  51  mm  thick  compact  'C’-shaped  specimens  (51CC)  and  the  51  mm 
thick  compact  tension  specimens  (51CT)  the  P^x/Pj  ratios  were  approximately 
1.12.  It  is  considered,  therefore,  that  the  Kg  values  reported  in  Table  2 
approach  the  Klc  of  the  material.  Thus  in  the  transverse  orientation  Klc  is 
approximately  130  MN  m-3  . 
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TABLE  1 


Chemical  Composition  of  steel  from  a  4.5  inch  Naval  Gun 

<wt%) 


c 

Mn 

P 

S 

Si 

Ni 

cr 

Mo 

0.32 

0.62 

0.009 

0.011 

0.20 

2.56 

0.84 

0.57 

TABLE  2 


Mechanical  Properties  of  Gun  Steel 


0.2 


GTS 


Kg  25CT  (transverse) 

Kg  5 ICC  (transverse) 

Kg  5 1CT  (longitudinal) 


79i'mn  m-2 


932  MN  m 


-2 


-3/2 


133  MN  m 
130  MN  m-3/2 

150  MN  m-3/2  (160  MN  m“3/2  est.  from  cod) 


Note:  All  properties  are  measured  at  20  c. 


Using  the  computer-interactive  test  method  described  above,  a  series 
of  standard,  single  specimen  J-integral  tests  conforming  to  astm  E  813  (2)  was 
undertaken  and  the  inferred  and  measured  crack  length  data  obtained  from  each 
test  are  given  in  Tables  3a  and  b.  Table  3a  summarizes  the  results  from  the 
actual  test  as  controlled  by  the  program  on  the  PDPll/04,  and  this  program 
incorporated  the  unloading  compliance  obtained  using  the  finite  element 
program,  PAFEC,  (17a, bl  (see  Figure  4).  Table  3b  summarizes  the  results  from 
a  re-run  of  the  experimental  data  using  JRERUN.BAS  on  the  VAX  l 1/780.  Errors 
in  the  measurement  of  aQ,  af  and  Aa  are  given  in  both  tables. 


TABLE  3a 


Summary  of  J- Integral  Test  Data  using  PAFEC-Generated  Compliance  Calibration. 
Data  from  PDPll/04  Control  Program,  1  (JCTS.BAS) 


specimen 

a 

O 

Err%<2) 

a 

f 

Aa 

NO. 

meas. 

mm 

CCBtpl. 

mm 

meas. 

mm 

comply 

mm 

Err% 

meas. 

rrrrt 

compl. 

mn 

Err% 

B1 

30.14 

30.74 

+2.0 

31.52 

31.77 

+0.88 

1.38 

1.03 

-25.4 

B2 

30.43 

31.07 

+2.1 

31.08 

31.61 

+1.7 

0.65 

0.54 

-16.9 

B3 

30.09 

30.53 

+1.5 

31.38 

31.76 

+1 .2 

1.29 

1.23 

-4.7 

B4 

30.18 

30.75 

+1.9 

30.47 

31.01 

+1.8 

0.29 

0.26 

-10.3 

B6 

30.50 

31.00 

+1.6 

31.13 

31.58 

+1.4 

0.63 

0.58 

.  -7.9 

B8 

29.94 

30.92 

+3.3 

- 

31.04 

- 

- 

0.12 

- 

BB13’ 

30. 18 

31.18 

+3.4 

- 

30.60 

- 

- 

0.0 

- 

B10(3) 

30.19 

30.69 

+1.7 

30.78 

31.18 

+1.3 

0.59 

0.49 

-16.9 

Bll 

30.27 

30.47 

+0.7 

30.40 

30.60 

+0.7 

0.13 

0.13 

+0.0 

Vote: 

(1) 

(2) 


Young's  modulus  -  203  GPa,  used  in  preliminary  calculations. 


Err%  is  defined  here  as  the  error  in  the  value  Inferred  from  the 
unloading  compliance  calibration. 


X  -  X 

Inferred  measured 


x  100 


measured 


Xmeasured  ls  ta)cen  as  absolutely  correct.  It  is  only  subject  to  errors 
in  physical  measurement  and  this  can  be  demonstrated  to  be  ±  .01  mm 
(±  10  *m) . 


No  heat  tinting  results  available  for  specimens  B8  and  B9. 
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TABLE  3b 


summary  of  J-lntegral  Test  Data  using  an  Experimental  compliance  Calibration. 
Data  from  VAXll/780  computer  Program,  ( JRERUN . BAS ) 


Specimen 

No. 

a 

meas. 

mm 

0 

compl . 

run 

Err%(l> 

a 

meas. 

rnn 

f 

compl. 

Btn 

Err% 

aa 

meas.  compl. 

ntn  mm 

Err% 

Bl 

30.14 

29.76(J*! 

wm 

31.52 

30.94{i5 

-1.8 

1.38 

1.18 

-14.5 

30.07  ^ 

Kfil 

31.22 

-1.0 

1.15 

-16.7 

B2 

30.43 

30.09 

m 

31.08 

30.68 

-1.3 

0.65 

0.60 

-7.7 

30.38 

WSm 

30.97 

-0.4 

0.59 

-9.2 

B3 

30.09 

29.54 

-1.8 

31.38 

30.77 

-1.9 

1.29 

1.23 

-4.7 

29.86 

-0.8 

31.06 

-1.0 

1.20 

-7.0 

B4 

30.18 

29.79 

-1.3 

30.47 

30.07 

-1.3 

0.29 

0.28 

-3.4 

30.10 

-0.3 

,  30.37 

k 

*-0.3 

0.27 

-6.9 

B6 

30.50 

30.04 

-1.5 

31.13 

30.65 

-1.5 

0.63 

0.61 

-3.2 

30.35 

-0.5 

30.94 

-0.6 

0.59 

-6.3 

B8<5) 

29.94 

29.96 

+0.1 

_ 

30.08 

_ 

_ 

0.12 

30.27 

+1.1 

30.39 

- 

0.12 

- 

09(5) 

30.16 

29.72 

-1.5 

_ 

29.72(®| 

_ 

0.0?(6) 

0.00 

__ 

30.04 

0.4 

30.04  ® 

- 

0.00 

- 

BIO 

30.19 

29.76 

-1.4 

30.78 

30.27 

-1.7 

0.59 

0.51 

-13.6 

30.07 

-0.0 

30.57 

-0.7 

0.50 

-15.3 

Bll 

30.27 

29.51 

-2.5 

30.40 

29.64 

-2.5 

0.13 

0.13 

1 

29.81 

-1.5 

29.94 

-1.5 

0.13 

1 

(1)  Err%  as  defined  for  Table  3a. 


(2) 

(3) 

(4) 


(5) 

(6) 


Upper  line  are  results  calculated  using  Young's  modulus  -  203  GPa. 
Lower  line  are  results  calculated  using  Young's  modulus  -  210  GPa 


This  value  i3  not  the  same  as  that  given  in  Table  3  of  [31,  because  in 
the  present  analysis  the  "next-to-f inal"  unloading  data  pairs  were  used 
for  linear  regression  analysis.  More  unloading  data  pairs  are  available 
and  hence  the  estimate  of  aa  is  possibly  more  accurate  than  in  [31 
despite  the  fact  that  a  small  amount  of  crack  growth  has  occurred 
between  unloads. 


Mo  heat  tinting  results  available  for  specimens  B8  and  B9. 
B9  showed  no  crack  growth. 
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V.s 


For  specimens  Bl  to  Bll  different  final  crack  extension  values  were 
attained  for  each  individual  test  which* means ,  therefore,  that  the  series 
could  also  be  treated  as  one  multiple  specimen  J  test.  The  final  craclc  length 
for  each  of  the  tests  was  marleed  by  heat  tinting;  some  of  the  tinted 
specimens  are  shown  in  Figure  11.  The  final  Aa  values  measured  from  the  heat 
tinted  specimens  are  included  in  the  summary  of  the  test  results  given  in 
.Tables  3a  and  b  and  were  measured  using  an  unweighted  ll  point  average*  which 
included  the  two  outer-surface  locations.  In  some  cases,  where  the  final 
craclc  extensions  are  quite  small,  less  than  0.20  mm,  it  is  possible  that  the 
requirements  for  a  valid  (or  even  useful)  J- integral  test  are  not  met  (see  (21 
section  9.3) . 

•• 

From  Table  3a  it  can  be  seen  that  the  original  PAFEC-generated 
compliance  calibration  used  in  the  PDPll/04  control  program  gives  what  appears 
to  be  good  estimates  of  aQ  and  af,  with  maximum  errors  of  up  to  +3.4%. 

Further,  the  PAFEC  compliance  calibration  results  in  consistent  over-estimates 
of  actual  crack-length;  thus  the  relative  errors  between  any  sets  of  aQ  and 
af  values  are  apparently  very  small.  However  the  calculations  of  Aa  (which 
use  aQ  and  af  values)  are  often  subject  to  very  large  errors;  indeed  many  of 
the  estimates  of  craclc  extension  using  the  present  PAFEC  compliance 
calibration  do  not  meet  the  requirements  of  the  ASTM  standard  for  the  J- 
integral  test  (21  where  "the  final  cral'ck  extension  value,  as  predicted  by 
single  specimen  techniques,  must  agree  with  the  averaged  heat  tint  value 
within  15% ".  There  are  two  possible  reasons  for  this;  (1)  that  the 
compliance  calibration  used  initially  in  the  control  program  on  the  PDPll/04 
is  not  entirely  appropriate  for  the  test  and  (2)  the  number  of  points 
discarded  from  the  start  of  the  unloading  for  the  linear  regression  analysis 
is  not  sufficient  to  avoid  all  non-linear  behaviour.  This  program  originally 
discarded  2  unloading  data  pairs.  These  data  can  be  re-analysed,  however, 
using  alternative  compliance  calibrations  and  varying  the  number  of  discarded 
data  pairs. 

Using  the  program  JRERUN.BAS  on  the  VAX1 1/780  (or  PDPll/04)  an 
alternative  compliance  calibration  is  tested.  The  only  other  compliance 
calibration  presently  available  for  the  MRL  J03  specimen  is  the  experimental 
unloading  compliance  calibration  generated  using  13  mm  thick  MRL  J03  specimens 
1161,  see  Figure  4.  Using  this  alternative  calibration,  and  discarding  4  data 
points  from  the  3tart  of  the  unloading,  gives  slightly  different  estimates  of 
aQ  and  af;  the  results  of  the  calculations  are  summarized  in  Table  3b. 

The  experimental  compliance  (161  tends  to  slightly  under-estimate 
the  values  of  craclc  lengths  and  this  can  be  altered  further  by  using  a  higher 
value  of  Young's  modulus,  although  the  value  for  the  gun  steel  used  in  the 
present  work  has  not  been  measured.  A  Young's  modulus  value  of  203  GPa  was 


It  should  be  noted  that  ASTM  E  813  uses  nine  inner  points  plus  a  weighted 
two  outer  points  for  the  calculation  of  mean  crack  length.  The  present 
unweighted  ll  point  average  appears  to  give  better  correspondence  with 
crack  length  inferred  from  compliance  than  the  9+2  point  average  used 
in  E  813. 


used  for  the  normalization  of  the  data  in  the  experimental  determination  of 
the  unloading  compliance  calibration  but  a  value  of  210  GPa  may  be  more 
appropriate  for  the  thicle  specimens  used  in  the  J-integral  testing  of 
specimens  Bl-Bll.  Both  Young's  modulus  values  are  used  for  the  calculations 
presented  in  Table  3b.  While  an  E  value  of  210  GPa  appears  to  give  a  closer 
correspondence  with  the  heat  tinting  data  there  is  a  greater  relative  error 
between  aQ  and  af  values  using  this  higher  E  and,  as  a  result,  slightly  larger 
errors  in  the  estimations  of  the  da  values.  Over-all  the  use  of  the 
experimental  unloading  compliance  data  tl6l  in  the  J-integral  test  offers  a 
better  correspondence  (smaller  Err%  values)  between  the  inferred  craclc  length 
values  and  those  measured  from  heat  tinted  specimens  than  is  achieved  using 
the  original  PAFEC-generated  unloading  compliance  calibration. 

Typical  J-da  plots  for  the  single  specimen  test  methods  are  shown  in 
Figures  12,  13  and  14.  These  plots  used  J-aa  data  pairs  generated  by 
JRERUN.BAS  on  the  VAXll/780. 

<. 

The  blunting  line  used  in  all  calculations  is  given  by  equation  16, 
and  is  included  on  all  J-aa  plots. 

/ 

As  mentioned  above,  the  series  of  specimens  Bl  to  Bll  can  also  be 
treated  as  one  multiple  specimen  J-integral  test  on  this  gun  barrel  steel. 

The  af  values  used  in  this  method  are  those  obtained  from  the  heat  tinting 
experiments  and  the  J  values  corresponding  to  each  final  da  value  were 
calculated  from  the  area  under  the  load/load-line-displacement  record  at  the 
start  of  the  final  unload  of  each  test.  The  J-aa  plot  for  these  data  i3  shown 
in  Figure  15. 

An  estimate  of  fracture  toughness  in  terms  of  stress-intensity 
factor,  K,  is  made  using  the  relation 


K  -  VE'  J_  (17) 

J  ic 

c 


where  Kj  »  the  fracture  toughness  derived  from  a  multiple  specimen  or 
c  single  specimen  J-integral  fracture  toughness  test. 

Jlc  -  value  of  J  at  the  Intersection  of  the  blunting  line  and  the 
linear  regression  through  admissible  J-aa  data  points. 

E'  -  the  value  of  Young’s  modulus  used  in  unloading  compliance 
calculations.  In  the  cases  of  relatively  thlclc  specimens 
this  E'  value  would  approach  E/(l-r  ),  where  r  is  Poisson's 
Ratio . 


The  following  table  summarizes  the  results  of  the  J-integral 
calculations  using  both  multiple  and  single  specimen  test  methods. 


TABLE  4 


m 


Summary  of  J-Integral  Test  Results  for  a  NICrMo 
Gun  Steel  at  20°c 


K 


Specimen  No. 


KJC(1) 

aa  (l) 

max 

JIc(2) 

KJc(2) 

aa  (2) 
max 

aa  (measured) 
max 

MMrf3/2 

urn 

kum-1 

-  MMti~3/2 

rrm 

rrm 

127.4 

139.6 

131.3 
95.6 

127.4 

125.6 


73 
78 
80 
84 
»'  76 


121.6 

126.2 

127.5 

130.4 

124.2 

120.3 


103.3 


121.5  ±  15.9 


122.6  ±  7.3 


Multiple  Specimen 


Jlc  -  93  KN/rn 


Kjc  -  137.4  MNtn 


Votes : 

(1) 


From  the  original  PDPll/04  data,  see  [31.  2  discards  from  the  start  of 

an  unloading;  PAFEC-generated  unloading  compliance  calibration; 

E  -  203  GPa  and  a  MerXle-Corten  Coeff  of  2.265,  (too  high). 

From  JRERUN.BAS  on  the  VAX  11/780.  4  discards  from  the  start  of  an 

unloading;  experimentally  determined  unloading  compliance 
calibration;  E  -  203  GPa  and  a  Merkle-Corten  Coeff  of  2.216,  (correct 
value) . 

No  heat  tinting  results  for  these  specimens. 


All  data  for  aa 


(viz  B1,B2,B3,B4,B6,B9  ( Aa-0 ) ,  B10,Bll). 


B 


'j  *  .•  ^  * 


The  results  of  the  calculations  of  JIC'  given  in  Table  4,  for  two 
methods  of  analysis  show  that  the  contrcrl  program,  JCTS.BAS  and  the 
arithmetically  more  accurate  JRERUN.BAS  produce  similar  Jlc  values,  some 
tests  such  as  B9  and  Bio  give  measurably  lower  Jlc  values  than  the  remaining 
bulk  of  the  tests.  The  reason  for  this  is  not  apparent  because  the  same  test 
method  and  data  analysis  was  applied  to  these  specimens,  and  hence  it  must  be 
concluded  that  specimens  B9  and  Bio  are  representative  of  the  lower  toughness 
values  in  what  may  be  normal  scatter  for  this  material,  viz.  Kjc  -  122.6  ± 

7.3  MN  m-3. 

From  the  results  summarized  in  Table  4  it  can  be  seen  that  the 
multiple  specimen  test  gives  a  higher  Jlc  value  than  the  series  consisting  of 
single  specimen  tests.  This  difference  arises  because  of  distribution  of  the 
J-Aa  data  points  between  the  0.15  and  l.so  mm  exclusion  lines.  The  multiple 
specimen  test  has  these  points  evenly  distributed  between  the  exclusion  lines, 
but  in  the  cases  of  the  single-specimen  tests  these  data  are  often  closer  to 
the  blunting  line,  thus  producing  lower  Jlc  values.  Despite  limited  multiple 
specimen  J-integral  data  both  the  single  specimen  and  multiple  specimen  KJc 
values  agree  well  with  the  known  Kg  values  of  the  gun  steel.  Unfortunately, 
there  are  no  Klc  values,  but  because  the  P^/Pj  values  for  the  test  are 
approximately  1.12  (261,  then  the  comparison  between  the  Kjc  values  measured 
in  the  present  work  and  the  known  Kg  v&lues  is  not  unreasonable. 

These  preliminary  Jlc  tests  using  the  gun  steel  highlight  some 
deficiencies  in  the  computer-interactive  single  specimen  J-integral  test 
method,  but  these  can  be  overcome  by  the  use  of  JRERUN.BAS.  The  first 
deficiency  is  that  JCTS.BAS  should  discard  more  than  two  unloading  data  pairs 
for  the  linear-regression  analysis  since  there  is  evidence  of  non-linear 
behaviour  within  the  first  four  data  pairs.  This  contributes  to  errors  in 
crack  length  measurement,  although  this  would  not  readily  be  apparent  in  the 
"rounding  off"  of  data  to  2  decimal  points.  Another  deficiency  is  that  within 
JCTS.BAS  the  correlation  coefficient  for  the  linear  regression  analysis  should 
be  calculated  and  documented.*  At  the  present  stage  of  development  this  is 
only  handled  by  JRERUN.BAS.  a  further  limitation  of  the  JCTS.BAS  program  is 
its  lack  of  generality;  certain  lines  within  the  program  have  to  be  edited  if 
test  parameters  are  changed;  e.g.  discards  of  data  pairs  in  the  linear 
regression,  compliance  calibration  and  displacement  calibration.  At  a  later 
stage  of  development  the  program  could  be  expanded  for  these  to  be  written 
into  the  program  interactively. 


Those  factors  which  influence  the  accuracy  of  the  measurements  and 
calculations  are  discussed  in  Section  10. 


* 


Limited  TPA  of  PDP  ll/04  does  not  permit  these  calculations  to  be 
undertaken. 


9.  FURTHER  APPLICATIONS  OF  THE  COMPUTER-INTERACTIVE , 
SINGLE  SPECIMEN  J-INTEGRAt  TEST  METHOD 


At  its  present  level  of  development  this  test  method  has  been  used 
for  the  measurement  of  the  fracture  toughness,  Jlc,  of  6  mm  thick  MRL  J03 
‘specimens  of  a  precipitation  hardened  martensitic  stainless  steel,  STA60 , 

[211.  The  results  presented  in  C  2 7  J  were  found  to  behave  similarly  to  those 
data  reported  for  17/4  PH,  (221.  The  greatest  difficulty  encountered  with  the 
test  method  was  that  the  unloading  compliance  calibrations  used  in  the  work, 
PAFEC -generated  in  JCTS.BAS  and  experimentally  measured  in  JRERUN.BAS,  could 
not  accommodate  the  large  degree  of  crack  tunnelling  which  occurred  in  the 
testing  of  the  STA60.  This  problem  is  one  which  has  not  been  addressed  by  the 
developers  of  unloading  compliance  J-integral  test  methods. 

A  second  application  of  the  computer-interactive,  single  specimen  j- 
integral  test  method  was  in  the  measurement  of  the  fracture  toughness,  Jlc,  of 
a  105  mm  gun  tube  (non-autofrettaged)  of  NiCrMo  steel.  In  this  work  (281  it 
was  possible  to  determine  the  plane-strain  fracture  toughness,  KIc,  of  the 
steel,  with  which  to  compare  the  Jlc  value  obtained  from  JCTS.BAS,  as  KJ<;. 

The  agreement  between  Klc  and  Kjc  was  .particularly'  good.  Crack  tunnelling  in 
this  series  of  tests  was  minimal,  and  hence  agreement  between  inferred  and 
measured  Aa  values  generally  conformed  to  ASTM  E  813  (21. 

These  applications  of  the  computer-interactive,  single  specimen  J- 
integral  test  method  have  demonstrated  its  suitability  for  the  testing  of  high 
strength  steels. 


10.  DISCUSSION 


The  limited  TPA  of  the  PDPll/04  used  in  the  present  experiments 
meant  that  error  and  accuracy  checks  as  specified  in  (21  could  not  be 
undertaken  during  the  test.  Thus  most  of  the  accuracy  checks  have  to  be 
undertaken  retrospectively  using  the  JRERUN.BAS  program  and  this  makes  the  J- 
integral  test  method  developed  in  this  work  similar  to  the  KIc  test  1 1 ] ,  where 
a  completed  test  can  be  declared  invalid  on  the  basis  of  retrospective 
analysis. 

The  standard  (2)  requires  that  the  initial  linear-elastic  slope  be 
used  as  a  guide  for  the  accuracy  of  the  initial  test  set-up.  If  the  PDPll/04 
control  program  were  to  check  the  loading  compliance  following  the  elastic 
load-up  (Z9-1),  then  it  would  be  essential  for  the  elastic  load  limit  of  the 
test  to  be  set  to  the  maximum  final  pre-cracking  load  of  the  specimen.  This 
means  that  if  the  test  were  found  to  be  progressing  incorrectly  then  it  could 
be  stopped  and  re-started  without  the  need  for  further  pre-cracking,  i.e.  the 
specimen  would  still  satisfy  the  final  pre-cracking  requirements  for  a  valid 
te3t.  Unfortunately,  this  load  is  relatively  low.  For  example,  in  the  case 
of  the  gun  steel  (Section  8)  the  elastic  load  limit  would  be  25  kN  or  lower! 


In  the  present  method  the  specimens  are  checked  for  dimensional 
accuracy  and  so  it  can  be  assumed  that  *if  the  clip  gauge  seats  correctly  in 
the  specimen  then  errors  in  the  load/load-line-displacement  record  will  be 
minimal.  (Correct  seating  of  the  clip  gauge  can  be  checked  using  very  small 
pre-loads  and  unloads.)  The  alignment  of  the  test  machine  is  also  checked 
prior  to  starting  a  series  of  j-tests.  It  is  considered,  therefore,  that  the 
preliminary  check  of  the  compliance  is  not  necessary.  Any  discrepancies  in 
compliance  with  the  experimentally-determined  value  can  probably  be 
accommodated  by  small  changes  in  Young's  modulus  in  JRERtfN.BAS  (see  below). 
This  check  of  the  initial  linear  elastic  compliance  is  done  indirectly  by  the 
JRERUN.BAS  program.  It  can  be  shown  that  a  ±  7%  variation  in  normalized 
compliance  at  a/W  -  0.60  produces  a  ±  2.3%  variation  in  inferred  a/w  value, 
or,  for  a  25  mm  thick  specimen,  a  crack  length  of  30.60  mm  ±  0.70  mm.  with 
the  exception  of  specimens  B8  and  B9,  the  initial  crack  length  estimates  in 
Table  3(a)  fall  within  this  range  suggesting  that  most  tests  were  correctly 
set-up. 

The  area  calculation  check  in  the  linear-elastic  loading  region  is 
considered  more  important  than  the  linear-elastic  loading  compliance  (above) 
because  it  will  give  an  indication  of  the  errors  in  the  area  under  the 
load/displacement  curve  introduced  by  the  area  calculation.  This  can  be 
checked  further  against  the  analogue  plot  (see  below).  To  undertake  this 
numerical  check,  it  is  desirable  that  the  elastic  load  limit,  P7,  be  set  as 

high  as  possible  at  the  start  of  the  test  so  that  a  relatively  large  area  is 

measured  and  checked  against  a  calculated  area  (area  of  a  triangle  in  the 
elastic  region) .  in  the  present  form  of  the  test  control  program  this  elastic 

load  limit  is  usually  set  higher  than  the  maximum  fatigue  load  (for  the  final 

7.6  mm  of  crack  growth)  and  this  precludes  the  re-use  of  the  specimen  if  a 
test  does  not  proceed  for  any  reason.  Experience  with  steels  has  shown, 
however,  that  with  the  present  instrumentation  and  level  of  development  of  the 
test  method,  preliminary  compliance  checks  are  unnecessary,  and  that  once 
started  the  test  should  not  be  interrupted.  The  preliminary  J-lntegral  tests 
conducted  on  the  gun  steel  revealed  a  slight  error  in  the  calculation  of  the 
area  under  the  load/load-line-dlsplacement  record,  producing  an  over-estimate 
of  this  area  in  the  linear  elastic  region  up  to  load  P7  (see  Table  5 
below) .  This  has  been  modified  and  now  most  MRL  J03  specimens  give  an  error 
in  the  elastic  region  of  less  than  +2.0%  [27,281.  The  magnitude  of  error  is 
calculated  by  the  ACCESS. BAS  program  on  the  PDP  11/04.  This  particular  error 
check  gives  a  good  indication  of  the  accuracy  of  the  area  measurement  in  the 
initial  stages  of  the  test  where,  because  of  the  slope  of  the  load/load-line- 
displacement  record,  the  calculation  errors  should  be  the  largest  of  the  test. 
This  may  al30  give  a  measure  of  the  non-linearity  which  is  observed  at  the 
start  of  most  fracture  toughness  tests. 


A  check  against  the  area  under  the  analogue  load/load-line- 
displacement  records  for  specimens  Bl-Bll,  using  a  Kontron  Blldanalyse  Minl- 
MOP  for  area  determination,  shows  good  agreement  between  the  calculated  and 
measured  areas.  These  are  summarized  in  Table  5.  It  is  not  possible  to 
determine  errors  for  the  computer-calculated  areas  because  there  are  no 
measured  areas  which  are  absolutely  correct.  The  areas  measured  from  the 
analogue  plots  are  subject  to  undetermined  errors  (but  which  are  likely  to  be 
small  anyway).  These  are,  for  example,  the  accuracy  of  the  analogue  plot 
(pen  thickness,  linearity  of  scales  calibration,  etc.)  and  the  accuracy  of 
usage  of  the  Mini-MOP,  see  Table  5.  Within  these  limitations,  the  agreement 
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between  the  calculated  and  the  measured  areas  £or  specimens  Bl-Bll 
demonstrates  that  the  JCTS.BAS  program  produces  acceptable  estimates  of  area 
under  the  load/load-llne-dlsplacement  record  from  which  J  values  can  be 
determined . 


TABLE  5 


Area  Calculations  for  Single  Specimen  J-lntegral  Tests 


Specimen 

No. 

Measured  Area^1* 
kN  mm 

Calculated  Area 
kN  mm 

(2) 

Error  to  Elastic' 
Load  Limit 

% 

Elastic 

Load 

Limit 

kN 

B1 

35.97  ±  0.233 

34.92 

1.58 

35 

B2 

32.07  ±  0.035 

32.68 

r  2.3i 

35 

B3 

39.84  ±  0.093 

39*.  14 

4.60 

35 

B4 

23.20  ±  0.066 

23.63 

2.45 

35 

B6 

27.07  ±  0.072 

27.21 

2.05 

35 

B8 

18.42  ±  0.060 

19.51 

4.11 

N  ■ 

35 

B9 

13.90  ±  0.037 

19.25 

3.54 

30 

BIO 

29.71  ±  0.032 

29.97 

1.21 

35 

Bll 

19.85  ±  0.052 

18.82 

2.48 

35 

Votes : 

(l)  Measured  using  Mlni-MOP;  standard  deviations  reflect  reproducibility  of 
area  measurements  using  a  minimum  of  4  determinations. 


.  ,  measured  calc 

(2)  Err%  -  - r -  x  100;  results  from  ACCESS. BAS. 

Acalc 


One  limitation  of  the  test  procedure  at  the  present  level  of 
development  is  in  the  limited  number  of  load/load-line-dlsplacement  pairs  used 
for  the  calculation  of  the  unloading  compliance  (Z9-3).  Dependent  on  the 
displacement  calibration  of  the  system,  the  number  of  admissible  data  pairs 
should  be  in  excess  of  20  (see  (29}  for  a  discussion  of  unloading  compliance 
instrumentation) .  Some  computer-interactive  test  methods  use  more  than  60 
data  pairs  and  to  obtain  such  a  large  number  of  points  both  those  from  the 
unloading  and  the  re-loading  are  used  [ 3 o 1 .  At  the  present  time  it  is  not 
possible  to  determine  the  optimum  number  of  data  pairs  for  determination  of 


the  unloading  compliance  of  the  specimens;  however,  to  assist  in  the 
assessment  of  the  over-all  acceptability*  of  data,  the  linear  regression 
calculations  within  JRERUN.BAS  also  give  the  coefficient  of  determination,  r2, 
for  each  estimate  of  unloading  compliance,  ideally,  this  value  should  be  in 
excess  of  0.9999, [221,  however,  for  the  Bl-Bll  series  undertaken  in  this  work 
most  of  the  linear  regression  calculations  for  each  unloading  give  values  in 
excess  of  0.999.  A  typical  output  from  JRERUN.BAS  is  given  in  Appendix  5, 
where  it  can  be  seen  that  the  coefficients  of  determination  for  unloading 
compliance  data  are  in  excess  of  0.999.  If  an  increase  in  the  number  of  data 
pairs  for  the  unload  is  found  to  be  desirable  then,  in  the  cases  of  specimen 
geometries  for  which  the  present  J-integral  test  control  program  has  been 
developed,  this  can  be  achieved  simply  by  changing  the  span  control  of  the 
MTS-810  system  so  that  the  number  of  bits  from  the  MSG  results  in  a  smaller 
increment  in  stroke. 


»  '  -  « 
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Another  potential  source  of  error  in  the  single-specimen  J-integral 
test  method  is  in  the  non-linear  unloading  effects  which  occur  with  some 
materials.  This  can  introduce  significant  errors  in  the  calculation  of 
unloading  compliance  1291  and  also  small  errors  in  the  estimation  of  area 
under  the  load/load-line-dlsplacement  curve.  To  allow  for  these  non-linear 
effects  (time-dependent  recovery)  a  "hold"  in  the  displacement  is  often 
desirable  before  unloading  see  (22  and’  291.  The  material  tested  in  this 
preliminary  work  was  a  high  strength  steel  and  as  a  result  the  non-linear 
effects  at  the  start  of  an  unload  are  small.  This  can  be  seen  in  Figure  16, 
where  several  unloads  are  re-plotted  by  GRAPH  (24]  after  running  the  data 
through  ACCESS1.BAS  on  the  VAX  11/780.  It  is  noted  that  a  very  small  degree 
of  non-linearity  is  evident  at  the  top  (start)  of  the  unloading  record.  To 
overcome  the  possibility  that  any  data  pairs  which  may  constitute  the  non¬ 
linear  part  of  the  record  are  admitted  to  the  regression  analysis  the  control 
program  (JOTS  for  example)  now  discards  the  first  4  load/load-line- 
dlsplacement  data  pairs.  In  the  case  of  the  data  analysis  programs, 
JRERUN.BAS  on  the  VAX  11/780  and  the  PDP  11/04,  all  data  pairs  are  listed 
prior  to  linear  regression  analysis,  and  the  data  can  be  re-examined  and  the 
number  to  be  discarded  input  before  the  calculations  begin.  The  number  of 
data  pairs  which  are  discarded  can  vary  for  successive  unloads  in  the 
JRERUN.BAS  program  and  hence  is  listed  in  the  summary  of  the  calculations 
produced  by  the  program. 


Errors  in  the  estimations  of  Aa,  and  to  a  lesser  extent  the 
calculation  of  J,  can  also  arise  as  a  result  of  electrical  and  mechanical 
hysteresis  in  the  load  and  displacement  measurement  systems,  using  a  250  kN 
load  cell,  the  load  can  be  measured  to  within  ±  1.0%  and  the  displacement  clip 
gauge  has  a  maximum  error  of  ±  0.1%.  Electrical  noise  is  minimised  because 
each  load/load-line-displacement  data  pair  is  an  average  of  10  readings  taken 
over  0.01  second. 


Looking  back  to  the  data  in  Tables  3(a  and  b),  generated  by 


/ 


The  acceptability  of  the  data  can  be  tested  quantitatively  on  the  basis  of 
the  accuracy  of  predicted  versus  measured  aQ  and  af  values. 
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JRERUN .BAS,  it  can  be  seen  that  the  agreement  between  measured  and  calculated 
crack  length  values  is  closer  than  those  obtained  using  the  PAFEC-generated 
compliance  calibration  for  the  MRL  J03  specimens.  The  sensitivity  of  the 
calculated  results  to  the  compliance  calibration  and  Young's  modulus  value 
used  for  the  normalization  of  raw  compliance  data  has  been  demonstrated  in 
Section  8,  and  is  clearly  an  important  aspect  of  the  appraisal  of  the  test 
■procedure . 


It  is  difficult  to  assess  the  accuracies  to  which  the  J-integral 
value  and  corresponding  Aa  value  are  measured.  All  aspects  of  the 
calculations  and  error  analyses  must  be  documented  during  a  test  or  in  the  re¬ 
run  of  the  raw  data  through  JRERUN.BAS. 

The  final  point  to  be  discussed  is  the  accuracy  to  which  the  Jlc 
value  can  be  determined  from  a  set  of  admissible  J-Aa  data  pairs.  It  is 
usual  to  fit  a  straight  line  through  all  admissible  data  pairs  whereby  Jlc  is 
calculated  at  the  point  of  intersection  of  this  linear  regression  line  and  the 
blunting  line  [21.  This  method  has  been  used  in  the  present  work.  However, 
it  is  instructive  to  re-examine  some  of  the  results  given  in  Table  4  in  more 
detail.  The  Jlc  values  under  consideration  are  those  obtained  from 
JRERUN.BAS  using  the  experimentally-determined  unloading  compliance 
calibration.  The  calculations  for  the  Jlc  values  are  summarized  in  Table  6. 

From  the  ±  95%  confidence  limits  on  the  estimation  of  JIc  for  each 
of  the  specimens  tested  in  this  work  it  can  be  seen  that  3ome  of  the  Jlc 
values  are  poor  estimates  of  fracture  toughness.  Appreciable  errors  in  the 
estimation  of  Jlc  arise  when 

(i)  the  number  of  admissible  J-Aa  data  pairs  is  low. 

(li)  the  J-Aa  data  pairs  are  grouped  away  from  the  0.15  mm 

exclusion  line  and  towards  the  1.50  mm  exclusion  line,  [21. 

(ill)  the  J-Aa  data  pairs  exhibit  scatter  within  a  single  specimen  test. 


It  is  clear  that  the  ideal  grouping  of  the  J-Aa  data  points  for  the 
purposes  of  the  linear  regression  analysis  is  equi-spaced  between  the 
exclusion  lines  and  that  the  number  of  data  pairs  should  be  in  excess  of  12. 
This  is  not  satisfied  by  some  of  the  tests  conducted  in  this  work.  where  the 
above  criteria  are  met,  as  in  the  cases  of  specimens  Bl  and  B2,  the  ±  95% 
confidence  interval  is  small.  The  data  for  specimen  B4  have  been  analysed  by 
successively  using  more  J-Aa  data  in  the  regression  analysis.  unfortunately, 
where  6  and  9  data  points  are  used  the  extra  data  points  are  to  the  left  of 
the  0.15  mm  exclusion  line.  However  the  result  demonstrates  the  need  for 
sufficient  data  points  of  optimum  spacing.  Whether  or  not  the  scatter  of 
J-Aa  points  can  be  reduced  for  a  single  specimen  J-integral  test  is  difficult 
to  assess.  If  all  physical  aspects  of  the  J-integral  test  are  correct  (see 
above)  then  any  scatter  observed  in  the  test  data  must  be  regarded  as  a  result 
of  the  response  of  the  material  to  the  loading  (and  unloading)  of  the  J- 
integral  test  procedure.  The  effects  of  scatter  on  the  estimation  of  Jlc  can 
be  reduced  by  keeping  the  number  of  admissible  data  pairs  high  and  this  can  be 
done  quite  simply  using  computer-interactive  test  methods. 


TABLE  6 


Summary  of  Linear  Regression  Calculations  for  Jlc  Determinations 


Specimen 

No. 

Jlc  value5 

JcN  m-1 

No.  of 
Admissible 
Data  Points 

Y  Intercept 
JcN  m-1 

±  95%  Confidence  2 
limits 

JcN  m-1 

Bl 

72.8 

13 

70.2 

±  5.5 

B2 

•  78.4 

9 

73.0 

±  9.1 

B3 

80.1 

ll 

71.0 

±  22.7 

B4 

not  calculated 

5 

41.9 

±  75.3 

B4 

83.8 

6 

79.5 

±  32.2 

B4 

not  calculated 

9 

48.8 

z' 

±  8.8 

B6 

75.9 

72.6 

±  6.8 

B8 

71.3 

no  values 

- 

- 

B9 

>  60 

no  values 

- 

- 

BIO 

57.4 

9 

55.4 

±  22.5 

Bl  1 

>  70 

no  values 

- 

- 

(1)  On  the  blunting  line 

(2)  On  the  Y  axis 


In  summary  it  can  be  seen  that  most  of  the  errors  discussed  above 
are  approximately  the  same  for  each  test  using  the  MRL  J03  specimen.  The 
calculated  error  in  the  area  under  the  elastic  portion  of  the  load/load-line- 
displacement  curve  is  consistently  <+2%  (but  never  negative),  aQ  and  af  are 
usually  estimated  to  within  ±  2%  (but  if  similar  unloading  compliance 
calibrations  are  used  the  test-to-test  variations  in  craclc-length  estimates 
are  even  smaller)  and  Aa  is  usually  estimated  to  within  ±  15%.  Thus  under 
most  circumstances  the  computer-interactive  single  specimen  J-integral  test 
U3ing  MRL  J03  specimens  satisfy  the  standard  test  method  for  the  measurement 


11.  CONCLUSIONS 


A  computer-interactive,  single  specimen  J-integral  test  method  has 
been  developed  for  the  PDPll/04-MTS  810  system,  and  tested  using  the  mrl  J00 
specimen.  The  material  under  test  was  a  Ni-Cr-Mo  4.5  inch  gun  steel. 


(1)  The  control  program  on  the  PDP11/04  has  been  shown  to  conduct  a 
single  specimen  J-lntegral  test  which  produces  Kjc  values  which 
agree  well  with  lenown  plane-strain  fracture  toughness  values. 

* 

(2)  The  inferred  increments  in  cracle  length,  Aa,  obtained  using  the 
unloading  compliance  method  incorporated  in  the  control  program  are 
usually  found  to  agree  within  +15%  of  that  obtained  by  direct 
measurement.  Alternative  calibrations  and  slight  changes  to  the 
Young's  modulus  value  used  in  the  calculations  can  be  used  to 
improve  the  correspondence  between  measured  and  inferred  values  of 
cracK  length. 

(3)  The  computer-interactive  single  specimeiy'J-integral  test  conforms  to 
the  requirements  of  ASTM  standard  E  813  for  the  measurement  of  JIc. 
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FIGURE  1 


Schematic  representation  of  crack-tip  stress  and  strain  fields 
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FIGURE  2 


ASTM  E  *13  method  for  the  measurement  of  J 
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FIGURE  4 


An  experimental  unloading  compliance  calibration  for  an  MR1>  J03 
specimen,  data  from  [16]. 
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FIGURE  5 


An  analytical  unloading  compliance  calibration  for  an  MR L  J03 
specimen,  data  from  [17a] . 


FIGURE  7 


Flow  diagram  for  a  computer  interactive,  single  specimen  J- 
integral  test. 
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Loading,  unload  and  re-load  ramps  in  a  single  specimen  J- 
integral  test. 


FIGURE  8 


FORCE.  kN. 


FIGURE  9 


f°rc*/lo*d-line-<ii8placement  analogue  output  for  a 
•ingle  specimen  J-integral  test. 
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FIGURE  10  Micrograph  of  a  transverse  section  through  a  4.5  inch  Naval  gun. 
1000X  21  nital  etch. 
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FIGURE  11  Heat  tinted  specimens  to  mark  final  crack  lengths  in  specimens 

from  a  4.5  inch  gun  barrel  of  NiCrMo  steel. 
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Data  from  JRERUN.BAS,  see  AppendixEZ 
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FIGURE  12 


A  single  specimen  J-Aa  plot  for  specimen  B1 
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FIGURE  13 


A  single  specimen  J-&a  plot  for  specimen  B3< 
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FIGURE  14 


A  single  specimen  J-Aa  plot  for  specimen  BIO 
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FIGURE  15 


The  multiple  specimen  J-Aa  plot  for  all  specimens,  B1-B11 
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APPENDIX  1 


COMPUTER  PROGRAM  LISTINGS  FOR  SINGLE  SPECIMEN 
J  INTEGRAL  FRAC TURE  TOUGHNESS  TESTS 


(a)  Program  listings  for  PDP  ll/04  computer. 

(l)  JCTS.BAS 

100  REM  .  .  .AUSTRALIAN  GOVERNMENT  DEFENCE  SCIENTIFIC  AND  TECHNOLOGY  ORGANIZATION 
110  REM  ...MATERIALS  RESEARCH  LABORATORIES 
120  REM  JCTS5.BAS 
130  TIME(10) 

140  EEXjMP  /  BUTNCS'.  LINE  2260  ) 

150  ERASE 
160  PRINT 

170  PRINT  'THIS  PROGRAM  CONDUCTS  A  SINGLE  SPECIMEN  J  TEST'  r 
180  PRINT  /  PRINT  'THE  TEST  COMPLIES  WITH  ASIM»  813  ’ 

190  PRINT  /  PRINT  ’NOTE:  READ  THE  QUESTIONS  CAREFULLY’ 

200  PRINT  ’  FOR  A  LISTING  OF  ALL  PARAMETERS  RUN  'JLIST.BAS* ' 

210  PRINT  ’<CR>  /  INPUT  ZS 

220  ERASE  /  PRINT  '•••••  MTS  ANALOG  SET  UP  ROUTINE  •••••'  /  PRINT 
230  PRINT  'ENSURE:-' 

240  PRINT  STRAIN  CHANNEL  IS  OOTWECTED  TO  COMPUTER’ 

250  PRINT  STROKE  CTRL  M,M  RANGE  MJST  BE  *-10  M4  . 

260  PRINT  SPECIMEN  OUT’ 

270  N7-N7M 

280  PRINT  HYDRAULICS  ON’ 

290  PRINT  PROGRAM  CTRL  IS  LOCAL' 

300  PRINT  /  PRINT  'MASTER  CLEAR  (QNL  4  MJST  -  000),<CR>’; 

310  INPUT  ZS  /  PRINT 

320  PRINT  'SET  PROGRAM  CTRL  TO  REACTE,  <CR>  '; 

330  INPUT  ZS  /  PRINT 

340  PRINT  'ADJUST  SET  POINT  TO  GIVE  ZERO  DISP  (CHJL3-000)’ 

350  PRINT  ••••••  WARNING  . 

360  PRINT  ’  RAM  WILL  MOVE  SO  ENSURE  CLEARANCE  WITH  STROKE  ZERO' 

370  PRINT  ’<CR>  ' ;  /  INPUT  ZS  /  PRINT 

380  PRINT  ’  INSERT  SPECIMEN  USING  STROKE  ZERO,  <CR> 

390  INPUT  ZS  /  PRINT 

400  PRINT  'SPAN  #1  TO  1000,<CR>  ';  /  INPUT  ZS  /  PRINT 
410  PRINT  'INSTALL  CLIP  GAUGE  (632.02C.20)  AND  ZERO  OUTPUT’ 

420  PRINT  'COARSE  ZERO  IS  CN  THE  INPUT  BOX  AND  FINE  ZERO  IS  ON  THE  AMPLIFIER  WIT' 

430  PRINT  ’<CR>  /  INPUT  ZS  /  PRINT 

440  PRINT  'SET  UP  X-Y  RECORDER  FOR  LOAD  /  DISP  RECORD’ 

450  PRINT  ’USING  CALIBRATED  10V  POT  ON  442  CONTROLLER’ 

460  PRINT  ’  1AM  -  25  CM  (HORIZ)  -  5.0  V 
470  PRINT  ’<CR>  VMEN  SET  UP’  /  INPUT  ZS 

480  ERASE  /  PRINT  ’•••••  RANGE  SELECTION  * . /  PRINT 

490  PRINT  'LOAD  RANGE  (KN)’;  /  INPUT  RO 

500  R0-2047/R0  /  REM  RO  NCW  -  BITS/KN 

510  PRINT  'ENSURE  75%  STRAIN  RANGE  IS  SELECTED, <CR>  ’; 


520  INPUT  ZS 

530  Rl-1011  /  REM  1011  BITS  IS  1.0VM  CN  THE  75%  STRAIN  RANGE 
540  FI *80  /  REM  RAMP  RATE-80 
550  QUIT 
560  ENABL 

570  ERASE  /  PRINT  'TO  STOP  TEST  <BREAK>  'S’  WILL  DUMP  THRU  LINE1 5460  ' 

580  PRINT  ’<CR>  TO  OCNTINUE’  /  INPUT  ZS 

590  ERASE  /  PRINT . *  n/A  ZERO  READINGS  . /  PRINT 

600  PO-O  /  DO-O 

610  GOSUB  1760  /  REM  DACQ 

620  PRINT  'AYE  LOAD  -  ' ;P0/( 10*R0 ) ; *  KN  ( OR  * ;P0/l0  ; '  BITS)’ 

630  PRINT  ’AVE  LPD  -  ’ ;DO/(lO*Rl ) ; ’  M  (OR  ';DO/lO;'  BITS)' 

640  PRINT  'ZERO  RESPECTIVE  AMPLIFIER.  REPEAT  READINGS  (Y/n)'; 

650  INPUT  ZS  /  IF  ZS-’Y'  THEN  600 

660  ERASE  /  PRINT  .  DATA  STORAGE  /  priot 

670  PRINT  /  PRINT  'ENSURE  DISK  IS  IN  DY1 :  ’ 

680  PRINT  ’SPECIMEN  TITLE’;  /  INPUT  TS 

690  OPEN  'DY1 :  ’♦TS*  ’  .DAT'  FOR  OUTPUT  AS  FILE  VF1  (20000) 

700  FOR  Z-0  TO  20000  /  VFl(Z)-0  /  NEXT  Z 

710  OPEN  'DY1 :  ’*TS* '  .PLT'  FOR  OUTPUT  AS  FILE  VF2(300) 

720  FOR  Z-0  TO  300  /  VF2(Z)-0  /  NEXT  Z  r 

730  OPEN  'LS:  '  PCX  OUTPUT  AS  FILE  #1 

740  ERASE  /  PRINT  . . ;TS ; '  DIMENSIONS  •••••*  /  print 

750  PRINT  'SPEC  WIDTH  (METRES)  -';  /  INPUT  W 
760  PRINT  ’SPEC  THICKNESS  (METRES)  -’;  /  INPUT  B 
770  PRINT  'E  (GPA)  -  * ;  /  INPUT  M 
780  PRINT  'PROOF  STRESS  (MPA)  -  ' ;  /  INPUT  U 
790  PRINT  'MERKLE-CORTEN  OO-EFFIC.  ' ;  /  INPUT  Q 
800  M-M*l .  00000E-09  /  REM  IN  N/M  2 

810  PRINT  /  PRINT  'TO  RETRIEVE  ’;TS;’  J  DATA, RUN  <ACCESS>’ 

820  COPY 

830  YF1  (0  )-B  /  VF1(1)-W  /  VFl(2)-M  /  VFl(3)-R0  /  VFl(4)-Rl 

840  PRINT  'ELASTIC  LIMIT  (KN)’;  /  INPUT  P7 

850  PRINT  ’MAX.  *  OF  UNLOADS';  /  INPUT  U7 

860  PRINT  ’#  OF  BITS  STROKE  INCR'  /  INPUT  M7 

870  VF2(l)-P7 

880  P7-INT(P7*R0 ) 

890  N-0  /  Nl-0  /  Z9-1  /  Z-0  /  D2-0  /  N7-0 

900  ERASE  /  PRINT  '**  RAMP  TO  ';P7/R0;'  KN  THEN  J  TEST’ 

910  PRINT  /  PRINT  ’X-Y  READY,  <CR>  TO  START  TEST’;  /  INPUT  ZS 
920  REM  PRINT  OF  DATA  FILE 

930  PRINT  #1  ,TAB(  10 )  ’PLCfTTING  DATA  FOR  THE  10%  UNLOAD  J  INTEGRAL  TEST  FOR  SPECIMEN  1  ;TS 
940  PRINT  #1  ,TAB( 5) ’READING  #’  ;TAB(20)  ’AREA’  ;TAB(35)  ’A/W’  ;TAB(50)  'DELTA 
CRKLNGIH'  ;TAB(75  )  ’  J  ’  ;TAB(85 )  'E5  ’  ;TAB(  100 )  'V' 

950  PRINT  #1  ,TAB(21 )  'KMA'  ;TAB(54  )  ’METRES  ’  ;TAB(72 )  'KN/M' 

960  REM  SET  UP  GRAPHICS 

970  ERASE  /  CXCN 

980  PHYL( 100,975,55,750) 

990  SCALE(0,0,2.00000E-0 3,0,350) 

1000  AXES(0 ,0 )  /  AXES (2. 00000E-03 , 350 ) 

1010  LABEL (TS* ’  DELTA  A,  METRES.’,’  J  VALUE ,KN/M’ , 5 . 00000E-04 ,25 .0  ) 

1020  INVEC  /  PLGT(O.O)  /  PLCT(200/(2*U*1000),250)  /  INVEC 

1030  INVEC  /  PLOT( 1 . 50000E-03 ,0 )  /  PLOT( ( 2 00/ ( 2 *U* 1 000 ) )♦ 1 . 50000E-03 , 250  )  /  INVEC 
1040  GOSUB  1130  /  REM  RAMP  ROUTINE 
1050  IF  Z9-3  THEN  GO  TO  1040 
1060  IF  Z9-5  THEN  GO  TO  1040 


1070  OOSUB  1870  /  REM  FIT  ROUTINE 
1080  MMtK('*’,YF2(5*Nl*3),J ) 

1090  PRINT 

*1  ,TAB(7);VF2(5*N1);TAB(22);VF2(5*N1M);TAB(37);VF2(5*N1»2);TAB(52);VF2(5,N1*3);TAB(70);VF2(5*N1- 
4 ) ;TAB(85 ) ;E5 ;TAB( 100 ) ;V 

1100  IF  VF2(5*N1*3)>-1 . 80000E-03  THEN  2260  /  REM  OUT  OF  PROG 
.1110  IF  N7>-U7  GO  TO  2260 
1120  GO  TO  1040 

1130  REM  •*  RAMP  SUB-ROUTINE  ** 

1140  IF  Z9-S  THEN  Z9-2 
1150  DIM  S (3 ) 

1160  CN  Z9  GO  TO  1180,1340,1500,1630  ■» 

1170  REM 

1180  REM  RAMP  STROKE  UP  CN  LOAD  WATCH  (Z9-1 ) 

1190  G3-So  /  REM  SO  IS  PRESENT  STROKE  POSITN 
1200  *n-0  /  DO-O 

1210  S-2  /  S(0)-0  S(1)-F1  /  S(2)-S9 
1220  FG(l ,S, 1 ,7,4) 

1230  PC*.  Z-0  TO  1000  /  NEXT  Z  /  REM  WAIT 
1240  N-N-l 

1250  GOSUB  1760  /  REM  EttOQ  r 

1260  VF1  (5*N)-N  /  VF1  (5*N*T  )-P5  /  VF1  (5*N*2  )-D5  /(,VF1  (5*NM)-Z9 
1270  IF  N-l  THEN  VF1  (5*N*3 )-P5  /  GO  TO  1300 
1280  P8“((P5+P4 )/ 2 ) 

1290  VFl(5*N+3)-VFl(5*N-2)*(P8)*((D5)-D7) 

1300  D7-D5  /  P4-P5 
1310  IF  P5>-P7  THEN  1740 
1320  S9-S9-1  /  GO  TO  1200 
1330  REM 

1340  REM  RAMP  STROKE  UP  ON  DELTA  DISPL.  WATCH  (Z9-2) 

1350  D9-D5+M7  /  REM  D9-DISPLACEMENT  END  LEVEL 
1360  PO-O  /  DO-O 

1370  S-2  /  S(0)-0  /  S ( 1 )-Fl  /  S(2 )-S9 
1380  FG(1,S,1,7,4) 

1390  FOR  Z-0  TO  1000  /  NEXT  Z  /  REM  WAIT 
1400  N-N*l 

1410  OOSUB  1760  /  REM  TO  DAGQ 
1420  P8-((PS*P4 )/ 2 ) 

1430  VF1  (5*N)-N  /  VF1(5*N*1)-P5  /  VFl(5*N*2)-D5  /  VF1  (5*N*3  )-VFl  (5*N-2)*(P8  )*(  (D5  )-D7 ) 

VF1 (5*N*4 )-Z9 
1440  A3-VF1 (5*N*3) 

1450  D7-D5  /  P4-P5 

1460  N5-N  /  REM  LAST  READING  NUMBER  AT  LOAD  UP 
1470  IF  D5>-D9  THEN  1740  /  REM  D9-NEW  DISP  LIMIT 
1480  S9-S9-1  /  GO  TO  1360 
1490  REM 

1500  REM  RAMP  DCMN  CN  LOAD  WATCH  (Z9-3 ) 

1510  N7-N7M 

1520  P9-.9*(P5)  /  D9-D5  /  REM  D9  IS  DISP  AT  START  OF  UNLOAD 
1530  PO-O  /  DO-O 

1540  S-2  /  S(0)-0  /  S(l )-Fl  /  S(2)-S9 
1550  FG( 1 ,S , 1 ,7 ,4 ) 

1560  FOR  Z-0  TO  1000  /  NEXT  Z  /  REM  WAIT 
1570  N-N*l 

1580  GOSUB  1760  /  REM  DACQ 

1590  VF1  (5*N)-N  /  VF1  (5*N*1  )-P5  /  VF1  (5»N-2  )-D5  /  VF1  (5*N*3)-A3  /  VF1  (5*N-4  )-Z9 
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1600 
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“•  • 

1630 
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1800 
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1840 

1850 
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* 
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1880 
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1900 

»  :  • 

1910 

r  *’ 

1920 

1930 

1940 

:  ''-J 

1950 

1960 

1970 

1980 

1990 

2000 

2010 

2020 

2030 

2040 

2050 

2060 

2070 

2080 

2090 

2100 

2110 

- 

2120 

2130 

2140 

2150 

IF  PJ<-P9  THEN  1740 
S9-S9-1  /  GO  TO  1530 
REM 

REM  RAMP  UP  TO  PRIOR  LEVEL  CN  DISP  WATCH  (Z9-4) 

PO-O  /  DO-O 

S-2  /  S(0)-0  /  S(1)-F1  /  S(2)-S9 
FG(l ,S,1 ,7,4) 

FOR  Z-0  TO  1000  /  NEXT  Z  /  REM  WAIT 
N-N*t 

OOSUB  1760  /  REM  DACQ 

VF1(5*N)-N  /  VFl(5*N*l)-P5  /  YF1  (5*N*2)-D5  /  VF1  (S*N*3  )-A3  /  VF1  (5*N*4 )-Z9 
IF  DS>-D9  THEN  1740  ' 

S9-S9*l  /  GO  TO  1640 
REM 

Z9-Z9-1  /  REM  RAMP  SWITCH 
RETURN 

rem  ••  d\cq  •• 

DIM  P(10).D(10) 

FOR  Z-0  TO  10  /  P(Z)-0  /  D(Z)-0  /  NEXT  Z 
DACQ(3,P,0,10)  /  DACQ(3  ,D,  1 , 10  ) 

START  ^ 

IF  P<10  THEN  1*10  ( 

IF  D<10  THEN  1*20 
QUIT 

FOR  Z-0  TO  9  /  P0-P0*P(Z)  /  D0-D0*D(Z)  /  NEXT  Z 

P5-PO/10  /  D5-D0/10 

RETURN 

REM  STRAIGHT  LINE  FIT 

X0-. 0748577  /  X1-. 017*943  /  X2— 2. 10639E-04 
X3-1 . 50063E-06  /  X4— 3.78693B-09 

K-R0/ (R1  *  1 . 00000E-06 )  /  REM  SCALING  FACTOR  FOR  VLL/P  IN  VM/KN 

REM  SET  SUMMATIONS  TO  ZERO 

El-0  /  E2-0  /  E3-0  /  E4-0  /  E5-0  /  E6-0  /  F-0 

IF  Nl>-1  THEN  2000 

RBI  USING  DATA  FROM  1ST  LOAD 

FOR  C-INT( . 9*N)  TO  INT(.l*N)  STEP  -1 

El -El *VF1 ( 5 *C* 1 )  /  E2-E2*VFl(5*C*2)  /  E3-E3-VF1 (S*C*l  )*VF1  ( 5*C*2  ) 

E4-E4-VF1 (5*C*2)  2 
F-F*l 

NEXT  C  /  GO  TO  2060 

REM  USING  DATA  FROM  10%  UNLOAD 

FOR  C-(N5*5)  TO  N  STEP  M 

El-El*VFl(5*C»l)  /  E2-E2*VFl ( 5*C*2 )  /  E3-E3*VF1 ( 5*C* 1 )*VF1 ( 5*C*2 ) 
E4-E4*VFl(5*C*2)  2 
F-FM 
NEXT  C 

REM  CALC  OF  SLOPE 

E5- (E3- ( (El *E2 )/F) )/(E4- (E2  2 )/F) 

V-  (K*M*B)/E5  /  REM  NORM.  OCMPL. 

R-X0*X1*VX2*V  2»X3 *V  3»X4»V  4 
Y1  -R*W 

IF  Nl>-2  THEN  2210 

VF2(3)-Y1 

IF  Nl>-I  THEN  2210 

VF2  (0  )-Yl  /  REM  ORIGINAL  CXAC3C  LENGTH 
REM  CALC  OF  Y  INTERCEPT 
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2160  E6“El/F-E5  *E2/F 

2170  RBrf  D8  IS  DISPL  AFTER  ALLOW  FOR  NONLINEAR  EFFECTS 
2180  D8—  E6/E5 

2190  REM  PRELIM  ESTIMATE  OF  AREA  (ELASTIC) 

2200  VF2(2)-(VFl(J*N*2)-D8)*VFl(S*NM)/2 
2210  N1  *N1  ♦  1 

2220  REM  CALC  OF  J  (KN/M) 

2230  J-Q*VF1  (5*N*3)/(B* (W-Yl ) *R0 *R1  •  1 000 ) 

2240  VF2(5*N1  )*N  /  VF2(  S*N1*1  )-VFl  ( 5*N*3 )/(R0*RI )  /  VF2 (5*NI*2 )-R  /  VF2(S*Nl*3)-Yl-VF2( j) 
VF2(S*N1*4)-J 
2230  REIURN 

2260  DUAP  /  MSWFG(1,2)  ' 

2270  COPY 

2280  ERASE  /  PRINT  'EM)  OF  TEST  _ PROGRAM  BY  T.Y.ROSE  AND  D.S. SAUNDERS,  1984' 

2290  CLOSE  /  STOP  /  END 


(2)  JRERUN.BAS 


100  REM  . .  .AUSTRALIAN  GOVERNMENT  DEFENCE  SCIENTIFIC  AND  RESEARCH  CRGINIZATiaN 
110  REM  ...MATERIALS  RESEARCH  LABORATORIES 
•  120  REM  JRERDN. . .  THIS  PROGRAM  WILL  RERUN  THE  JTEST  DATA 
130  ERASE 

140  PPTMT  «*»»■*'**  A  ******  *  A  'AA-A-A-A-A*  A~A"A"A"A~A~A  A  A  *****  AAA/tJfcit 

150  PRINT  "JRERUN,  A  PROGRAM  TO  RE-RUN  THE  J  TEST  USING  THE  DATA  FTLES" 

160  PRINT  "FROM  THE  PDP11  _ EDIT  LINES  1650-1740  TO  CHANGE" 

170  PRINT  "THE  COMPLIANCE  CALIBRATION  EtR  THE  MEASUREMENT  OF  CRACK" 

180  PRINT  "LENGTH  FROM  THE  UNLOADS  (Z9-3 
IN  THE  TEST  DATA  FILE)  " 

190  PRINT  "  *'*  * -A-*'*-***'* *■*'** A-A-A-A-A A-ArAr  A-A-A'A A' A  A  A-A-AA'  A  A  AAAArArArAAAAAArAAAr^  ***** 

200  PRINT  "<CR>  TO  CONTINUE  ’  INPUT  Z$ 

210  PRINT  "  A’AAArirAAAAAAAllrilrA-ArArA-*  *A*-AC**AA-AA'A»  AAA  AA-AArA-A'A-A'A'A-A *******  A  A'A'A  A  A  n 

220  PRINT  "SPECIMENS  FOR  WHICH  UNLOADING  COMPLIANCE  CALS.  ARE  AVAILABLE" 

230  PRINT  "SPEC  GEOMETRY"  ;TAB(  40 )  "TYPE" 

240  PRINT  "FRL  J03...PAFEC  1984";TAB<40)"1" 

250  PRINT  "FfiL  J03...EXPTAL  CCMP  ";TAB(40>"2"  /" 

270  PRINT  "  *****  *******  a  «  f 

1rk*****ic*r* ***** ********  A-A  A  A  A  A  A  « 

280  PRINT  "INPUT  SPECIMEN  TYPE  "  /  INPUT  Y 
290  CEf  Y  GOSUB  1650,1700 
300  PRINT  N$ 

310  PRINT  /  PRINT 

320  PRINT  'WHAT  IS  THE  OUTPUT  DATA  FILE  NAME  VIZ  NAME  '  /  INPUT  TS 
330  OPEN  'DY1:  '+T$+'  .DAT'  FOR  INPUT  AS  FILE  VFl(lOOOO) 

350  OPEN  'LS:'  FOR  OUTPUT  AS  FILE  #1 

360  PRINT  "HCW  MANY  READINGS  IN  THE  FTIE  "  /  INPUT  Q 

370  PRINT  "WHAT  IS  THE  MERKLE  CCRTEN  COEF  FOR  THIS  SPEC.  "  /  INPUT  M5 

380  DIM  P(40),L(40),A(40), 

N(40) 

390  DIM  R5(40),Q5(40),A5(40),U5(40),J5(40) 

400  DIM  F5(40),V5(40),C5(40),T5(40) 

410  Z9-2  /  N7-0 

420  PRINT  VF1  ( 3 ) ;  "BITS/KN” 

430  PRINT  VF1  ( 4 ) ;  "BITS/t*!" 

440  PRINT  "SPECIMEN  DIMENSIONS  SHOULD  BE  MEASURED  WITH  A  MICROMETER” 

450  PRINT  "SPECIMEN  THICKNESS  (METRES)  "  /  INPUT  B 
460  PRINT  "SPECIMEN  WIDTH  (METRES)  "  /  INPUT  W 
470  ERASE 

480  PRINT  "YOUNGS  MODULUS  IS  ”;VFl(2)  ;”PA” 

490  M-VFl(2) 

500  PRINT  "DO  YOU  WISH  TO  CHANGE  YOUNGS  MODULUS  Y/N  "  /  INPUT  S3 
510  IF  S$-"N"  GO  TO  540 

520  PRINT  "WHAT  IS  THE  NEW  VALUE, . .  .IN  GPA  "  /  INPUT  M 

530  M-M*1.00000E+09 

540  PRINT  /  PRINT 

550  PRINT  "CHECK  OF  VARIABLES" 

560  PRINT  "B-"  ;B,  "W-”  ;W,  "t*-”  ;M,  "MERKLE-CCRTEN  COEJF=»"  ;M5 
570  PRINT  "OK?  <Y/N>”  INPUT  S$ 

580  IF  SSO”Y”  THEN  GO  TO  440 


590  N3-VFl(3) 

600  N4-VF1U) 

610  K-N3/(N4*1.00000E*06) 

620  FOR  1-2  TO  Q*1 
630  Z-VF1 (5*1-1) 

640  IF  Z-l  GO  TO  680 

650  IF  Z-2  THEN  N-0 

660  IF  Z-(Z9-1)  THEN  OOSUB  880 

670  IF  Z-3  THEN  OOSUB  710 

680  NEXT  I 

690  IF  Z-3  THEN  GO  TO  880 

700  GO  TO  1290  ' 

710  REM  DATA  ACQUISITION 
FROM  POP  11  DATA  FILE 
720  IF  NOO  THEN  GO  TO  760 
730  PRINT 

740  PRINT  'UNLOADING  COMPLIANCE  DATA'  /  PRINT 
750  PRINT  'READING', 'LOAD', 'DISPL.'. 'AREA' 

760  N6-VFl( 5*1-5) 

770  P-VFl(5*I-4) 

780  L-VF1 (5*1-3)  r 

790  A-VF1 (5*1-2)  ( 

800  P(N)-P 
810  L(N)-L 
820  A(N)-A 
830  N(N)-N6 

840  PRINT  N(N),P/N3,L/to4,A/(N3*N4) 

850  N-N-l 
860  Z9-Z  /  N5-N 
870  RETURN 

880  PRINT  /  PRINT  /  PRINT 

890  PRINT  "INLQADING  COMPLIANCE  RERUN' 

900  REM  SET  SUMMATIONS  TO  ZERO 

910  El-0  /  E2-0  /  E3-0  /  E4-0  /  E5-0  /  E6-0 

920  E7-0  / 

E8-0  /  E9-0  /  Vl-0  /  F-0 

930  PRINT  "MAX.  NO  OF  UNLOADING  POINTS  IS  ';N5 
940  PRINT  TO*  MANY  TO  BE  DISCARDED  '  /  INPUT  0 

950  IF  (N5-0)<3  THEN  PRINT  'ERROR.  .  .NOT  ENOUGH  DATA  POINTS!"  /  GO  TO  940 
960  FOR  C-(0*O)  TO  (N5-l) 

970  E1-E1*P(C)  /  E2-E2*L(C)  /  E3-E3-P(C)*L(C) 

980  E4*E4-L(C)  2  /  E6-E6*P(C)  2 
990  F-F*l 
1000  NEXT  C 

1010  REM  CALC  OF  SLOPE 

1020  E8-E3-((E1*E2)/F)  /  E9-E4-((E2  2)/F) 

1030  E5-E8/E9 

1040  REM  CALC  OF  C0RR  COEFF 
1050  Vl-E6-((E1  2 )/F ) 

1060  E7-(E8  2)/(E9*Vl) 

1070  V-(K*M*B)/E5  REM  NO 

RM  COMPLIANCE 

1080  V5(N7)-V  /  C5(N7)-E7 

1090  R-X0  *X1 *V*X2  *Y  2*X3*V  3*X4*V  4 

1100  IF  N7>*1  GO  TO  1120 

1110  Y1-R*W  /  Y0-Y1  /  GO  TO  1130 
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1120  Y1-R*W 

1130  PRINT  "OffiCX  OF  UNLOADING  POINTS  "*;F;"  SHOULD  BE  ";N5-0 
1140  PRINT  "CRAOaENGTH  ;Y1 ;  ”VETRES” 

1150  PRINT  'READING' , "DATA  POINTS" , "NORM  CCMPL","C0RR  COEFF" , "AREA" , "aAT , "DELTA  A"."J  VALUE 
1160  J'Mi’ACcJ/CB’lW-Yl  )*N3*N4*1000 ) 

1170  U-N6-N5*l 

11*0  PRINT  U,F,V,E7,a(C)/(N3*N4 ) ,R,Y1  -Y0, J 
1190  Z9-Z 

1200  US(N7)-U  /  AS(N7)«A(C)/(N3*N4)  /  RS(N7)-R 
1210  Q5(N7)-Y1-Y0  /  J5(N7)-J  /  F5(N7)-F  /  T5(N7)-N5 
1220  PRINT  "<CR>  TO  CONTINUE  TO  NEXT  UNLOADING  SET"  /  INPUT  XJ 
1230  N7-N7*l  " 

1240  IF  I<Q  THEN  GO  TO  1280 

1250  PRINT  "LAST  UNLOADING  DATA  SET" 

1260  PRINT  *<CR>  TO  CONTINUE.  .PRINT  OUT  ON  LS : "  /  INPUT  XS 
1270  GO  TO  1290 
12*0  RETURN 
1290  CLOSE  YF1 

1300  PRINT  #1, ;TJ ;TS I"******" ;TS ;*••*•••" 

1310  PRINT  #1  /  PRINT  #1 

1320  PRINT  #1  .'SUMARY  OF  UNLOADING  COMPLIANCE  CALCS"/' 

1330  PRINT  #1,  "USING  PDPll/04  DATA  FILES  AND  RUN  CN  THE  PDPll/04" 

1340  PRINT  *1, 'FOR  FILE  ";TS 
1350  PRINT  #1.  /  PRINT  *1. 

1360  PRINT  *1  ."YOUNG 'S  MODULUS  FOR  THESE  CALCULATIONS  -",M;'PA" 

1370  PRINT  #1  ,'SPECIICN  THICKNESS  -";B;"kEIRES" 

13*0  PRINT  tl ,  "SPECIMEN  WIDTH  -"  ;W;"VCIRES" 

1390  PRINT  *1  /  PRINT  #1 
1400  PRINT  #1.NS 

1410  PRINT  *1.  "POLYNOMIAL  FIT  COEFFS" 

1420  PRINT  #1."X0-";X0;"  X1-';X1 X2'";X2 
1430  PRINT  ♦l."X3-";X3;"  X4-";X4;"  X5-";X5 
1440  PRINT  tl  /  PRINT  #1 
1450  PRINT  *1 . "VOKLE-CORTEN  CCEFF-  ";M5 
1460  PRINT  #1  /  PRINT  #1 

1470  PRINT  tl.  "READING  t  AT".  "TOTAL  NO.  OF",  "DATA  POINTS’ .  "NORM  CCMPL"  ,’OORR 
COEFF" ,  "AREA" .  "A/VT ,  "DELTA  A" ,  ’  J  VALUE" 

14*0  PRINT  tl  ."TOP  OF  IWLD',TAB(14)."DATA  POINTS"  ; TAB (2*  )  ."USED" 

1490  FOR  N-0  TO  N7-1 

1 500  PRINT  tl  ,U5  (N)  ,T5 (N)  ,F5  (N)  ,V5  (N)  ,C5  (N)  ,A5  (N)  ,R5  (N)  ,Q5  (N) .  J5  (n) 

1510  NEXT  N 

1520  PRINT  tl  /  PRINT  tl 
1640  GO  TO  1800 
1650  REM 

1660  NS-’J03  CAL.  FROM  PAFEC  19*4’ 

1670  X0-. 074*577  /  X1-. 017*943  /  X2— 2. 30639E-04 
16*0  X3*l .  50063E-06  /  X4  — 3 . 78693E-09 
1690  RETURN 
1700  REM 

1710  NJ-’J03  CAL.  FROM  CLARK  EXFT  DATA.  19*3  ’ 

1720  X0-. 106976  /  X1-. 0147547  /  X2  —  1 . 62423E-04 
1730  X3-8.75445E-07  /  X4— 1 . 77957E-09 
1740  RETURN 
1750  REM 

1*00  PRINT  tl.’OOOF  JRERUN.BAS .  .  .  PROOIAM  BY  D.  S.  SAUNDERS  AND  I. A. BURCH _ 19*6” 

1*05  CLOSE 
1810  EM) 


ai  .  a 


(3)  ACCESS. BAS 


100  RBI  . .  .AUSTRALIAN  GOVERNMENT  DEFENCE  SCIENTIFIC  AND  TECHNOLOGY  ORGINIZATICN 
110  REM  . .  MATERIALS  RESEARCH  LABORATORIES 

120  REM  ACCESS  PROC2AM  TO  ACCESS  DATA  IN  THE  HMM .DAT  FILE  AND  THE  #*##*#. GRA  FILE  FOR  THE  10% 
UNLOAD  i  TEST 

130  PRINT  WAT  IS  THE  DATA  FILE  NAVE?  *  INPUT  TS 
140  OPEN  *DV1 :  ’*TS* ’  .DAT'  FOR  INPUT  AS  FILE  VF1  (30000) 

150  OPEN  *LS: '  FOR  OUTPUT  AS  FILE  tl 

160  OPEN  'DY1:  ’♦TS*'.PLT'  FOR  INPUT  AS  FILE  VF2(500) 

170  PRINT  ’TEXT  REQUIRED?  Y/N'  INPUT  ZS  ' 

ISO  IF  ZS-'N'  GO  TO  270 

190  PR 

INT 

*1,TAB(5); . . . * . . . . . . . . . ...» 

200  PRINT  *1  .TAB  ( 5 ) ;  *  WWW  1 1  *THE  DATA  STCRFO  IN  THIS  PDP/ll  VIRTUAL  FILE  IS  IN  BINARY 

FORMAT******#*' 

210  PRINT  #1  ,TAB(5);’LQAD  THE  DISK  CN  TO  THE  VAX.  THEN  RLN  V.BAS  ON  THE  VAX.  THEN  THIS  WILL 
CREATE' 

220  PRINT  #1.TAB(5);'AN  OUT  PUT  FILE  OF  THE  ORIGINAL  PDP/ll  DATA,  AH)  IN  THE  SAME  LISTING 
FORMAT,"  f 

230  PRINT  *1  ,TAB(  5 ) ;  USE  THIS  OUTPUT  FILE  TO  CREATE  A  PLOTTING  DATA  FILE  TO  USE  WITH  THE  GRAPH  ’ 

240  PRINT  #1,TAB(5); 'PLOTTING  ROUTINE.  THE  DATA  CAN  BE  REPLOTTED  AM)  REANALYSED  USING  THIS 
ROUTINE' 

250  PRINT  *1  ,TAB(5);"MM»IHim»lll  THE  VAX  ROUTINES  ARE  l)  V.BAS  AH)  2)  ACCESS  1  .BAS 

mimuHiitHH- 

260  PRINT  *1 ,TAB(S) ;'*•*••*••••*••*•••*•**•*•••••*•••*•••••*•**•*•*•** . 


270  PRINT  *1  /  PRINT  *1 

280  PRINT  *1  ,TAB(20);TS;  ’  TEST  DETAILS' 

290  PRINT  *l.TAB(20)'SPBCn*N  THICKNESS  *  ;VF1  (0);  ’  M' 

300  PRINT  *1  .TAB (20 )  ’SPECIMEN  WIDTH’  ;VF1 ( 1 ) ;  ’  M' 

310  PRINT  #1 
320  R0-VF1(3)  / 

Rl-VFlU) 

330  PRINT  #1,TAB(20)  'SCALING  CONSTANTS  FOR  ’  ;TS ;  ’  TEST’ 

340  PRINT  *1  .TAB (25)  'FOR  LOAD,  R0-  ’;R0  ;’  BITS/KN’ 

350  PRINT  *1  ,TAB(25) 'FOR  DISPL.  Rl-  '  ;R1 ;  ’  BITSAM’ 

360  PRINT  *1 

370  PRINT  'DO  YOU  REQUIRE  A  PRINTOUT  FROM  THE  DATA  FILE?  (Y/N)  ’ 

380  INPUT  ZJ  /  IF  ZJ-’N’  THEN  740 

390  PRINT  'DO  YOU  WANT  TO  ACCESS  ALL  THE  DATA  Y/N’  /  INPUT  ZJ 
400  IF  ZJ-'Y'  GO  TO  440 

410  PRINT  "WAT  VALUE  DO  YOU  WANT  TO  START  AT?  '  /  INPUT  V 

420  PRINT  WAT  VALUE  DO  YOU  WANT  TO  FINISH  AT? '  /  INPUT  W 

430  00  TO  450 
440  V-I  /  W*2000 

450  PRINT  'THE  DATA  FILE  IS  STORED  IN  BITS’ 

460  PRINT  'DO  YOU  WANT  THE  DATA  FROM  '  ;TJ ; '  IN  KN.VW  AfC  KM4?  (Y/N)  ’  /  INPUT  YS 
470  IF  YS-'N’  THEN  610 
480  IF  XS-'Y’  THEN  345 

490  PRINT  #1  ,TAB(  18 )  'TEST  DATA  FOR  THE  10%  UNLOAD  J  INTEGRAL  TEST’ 

500  PRINT  #1  ,TAB(28 )  'DATA  FRCM  VIRTUAL  FILE  ’;TJ;’.EAT' 

510  PRINT  #1 


520  PRINT  #1,TAB(5) 'READING  *'  ;TAB(23)  'LOAD'  ;TAB(37)  'DISPL'  ;TAB(J0)  UNOORR.  AREA' ;TAB(6S ) 
SWITCH' 

530  PRINT  #1,TAB(2J)'KN’;TAB(J9)'VM';TAB(56)'WW' 

540  PRINT  *1 
550  FOR  Z-V 
TO  W 

560  IF  VFl(5*Z)-0  THEN  590 
570  PRINT 

*l.TAB(7);VFl(5,Z);TAB(2O);VFl(5*Z*l)/R0;TAB(35);VFl(5*Z*2)/Rl;TAB(S2);YFl(5,Z»3)/(R0,Rl) 

;VF1(5*Z»4) 

580  NEXT  Z 

590  PRINT  #1  /  PRINT  #1  /  PRINT  #1 
600  GO  TO  740 

610  REM  .THIS  PART  PRINTS  CUT  DATA  FRCN  THE  JTEST.DAT  VF 

620  PRINT  #1  ,TAB(  18 )  'TEST  DATA  FOR  THE  10%  UNLOAD  J  INTEGRAL  TEST’ 

630  PRINT  #1  ,TAB(32)  'DATA  FOR  '  ;TS ; '  TEST' 

640  PRINT  *1 . 

650  PRINT  *1,TAB(5) 'READING  *’ ;TAB(20 ) 'LOAD' ;TAB(35 ) 'DISPL' ;TAB(50 ) 'UNOORR. AREA’ ;TAB(65 )  ’ 
SWITCH' 

660  PR  IN 

T  *1  .TAB ( 20 ) 'B ITS '  ;TAB ( 36 ) ‘B ITS '  ;TAB ( 5 2 ) 'B ITS *B ITS '  ^ 

670  PRINT  #1 

680  FOR  Z-V  TO  W  ' 

690  IF  VF1 (5*Z)-0  THEN  740 
700  PRINT 

*1  .TAB(7);VFl(5*Z);TAB(22);VFl(S*Z*l);TAB(35);VFl(5,Z*2):TAB(52);VFl(5*Z*3);TAB(67);VFl(5 
710  NEXT  Z 

720  PRINT  *1  /  PRINT  #1  /  PRINT  #1  /  PRINT  *1 
730  GO  TO  840 

740  RB4  THIS  PART  PRINTS  OUT  DATA  FROM  THE  JTEST.GRA  VF 

750  PRINT  'THIS  PROGRAM  WILL  NOW  ONLY  PRINT  OUT  DATA  FROM  THE  J  CALCULATION  FILE  (VF2)  ’ 
760  PRINT  #1,TAB(20)' INITIAL  CRACK  LENGTH  FROM  ELASTIC  LOADING'  ;VF2(0 );  'M' 

770  PRINT  tl,TAB(20) 'ELASTIC  LOAD  LIMIT  FOR  THE  TEST  '  ;VF2(l ) ;  'RN' 

780  PRINT  *1  .TAB (20)  ’CRACK  LENGTH  CALCULATED  FROM  THE  FIRST  UNLOADING’  ;VF2 (3)  ;  'M' 

790  A1-VF2(6)  /  A2-VF2(2)/(VFl(3)*VFl(4)) 

800  F9-(A1-A2)*100/A2 

810  PRINT  #1,TAB(20) 'ERROR  IN  1ST  AREA  CALCULATION  AT  THE  ELASTIC  LOAD  LIMIT’ 

820  PRINT  #1  ,TAB(20) ’IS. . .  A(NCAS)-A(CALC)/a(CALC)  . . .  -  ’ ; F9 ;  ’ 

830  PRINT  tl  /  PRINT  *1 

840  PRINT  #1,TAB(  18)  'PLOTTING  DATA  FROM  THE  10%  UNLOAD  J  INTEGRAL  TEST' 

850  PRINT  #1  .TAB (32) 'DATA  FOR  ’;TS;’  TEST' 

860  PR 
INT  #1 

870  PRINT  *1,TAB( 5) 'READING  ♦’  ;TAB(20 )  ’AREA’  ;TAB(37  )  ’A/W’  ;TAB(53  )  'DELTA'  ;TAB(72  ) '  J  ' 

880  PRINT  01  ,TAB(52)  ’CRAQCLNGTH' 

890  PRINT  #l,TAB(20)'ra*4';TAB(53)’VeTRES’;TAB(7l)’ICN/M’ 

900  PRINT  #1 

910  FOR  Z-l  TO  100 

920  IF  VF2( J*Z)-0  THEN  940 

930  PRINT 

#I.TAB(7);VF2(S*Z);TAB(18);VF2(5*Z«1);TAB(36);VF2(5*Z*2);TAB(S0);YF2(5*Z*3);TAB(67);VF2(5 

940  NEXT  Z 

950  CLOSE  *1%  /  CLOSE  *2% 

960  PRINT  OI.'BC  OF  ACCESS.  BAS. ...  PROGRAM  BY  D.  S .  SAUNDERS  ’ 

970  EM) 


(b)  Program  listings  for  VAX  n/780  computer 


i 

f 

k 


(l)  V. BAS 


i 

i 


100  REM  . .  .AUSTRALIAN  GOVERNMENT  EEFQCE  SCIENTIFIC  AND  TECHNOLOGY  CRG3NIZATIQN 
110  REM  ...MATERIALS  RESEARCH  LABORATORIES 
*120  REM  PRCCSAM  TO  READ  BINARY  FILE  QF  THE  RT11  FORMAT 

130  INPUT  -INPUT  FTIENAME  [$$$9.DAT;#J  IE  A  PDP11  VIRTUAL  FILE  ",FH£NAMES 
140  INPUT  -OUTPUT  FUE  NAME  [5SS9.0UT;#]  :«  ",FHENAME29 
150  MAP  (A)  REAL  VT2<127) 

160  OPEN  FILENAME?  FOR  INPUT  AS  FILE  #1%,  £  _ 

VIRTUAL,  S 

RECCRDSIZE  512%,  £ 

MAP  A 

170  OPEN  FHENAME29  FOR  OUTPUT  AS  FILE  #2%,  RECCRDSIZE  132% 

180  K%-0% 

190  FOR  J%-1%  TO  79% 

200  FIND  #1%,  RECORD  J% 

210  GET  #1% 

220  FOR  1-0  TO  127  ^  . 

230  K%-K%+1% 

240  IF  K%05%  THEN  * 

PRINT  #2%,VT2(l) , 

ELSE 

K%-0% 

PRINT  #2%,VT2(I) 

END  IF 
250  NEXT  I 
260  NEXT  J% 

270  PRINT  "END  QF  V.BAS . PROGRAM  BY  D.WALKER  1984". 

280  END 
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( 2 )  JRERUN . BAS 


100  REM  . .  .AUSTRALIAN  GOVERNMENT  DEFENCE  SCIENTIFIC  AND  TECHNOLOGY  ORGANIZATION 
110  REM  ...MATERIALS  RESEARCH  LABORATORIES 

120  REM  JRERLN. .  .THIS  PROOtAM  WILL  RE-RIN  A  J  TEST  USING  TOE  DATA  FILES  FROM  TOE  PDP1 1 

130  PRINT  TAB(j ) ;"••*• . . . . . . . ••••••••• . ••••••" 

140  PRINT  TAB ( 3 ) ; ' JRERLN ,  A  PROOtAM  TO  RE-RLN  TOE  J  TEST  USING  TOE  DATA  FILES" 

150  PRINT  TAB ( 5 ) ;  'FROM  TOE  PDP11 . EDIT  LINES  1860-1960  TO  CHANCE" 

160  PRINT  TAB ( 5 ) ;  "THE  COMPLIANCE  CALIBRATION  FOR  TOE  MEASUREMENT  OF  CRACX" 

170  PRINT  TAB(5); "LENGTH  FROM  TOE  LNLQADS  (Z9-3  IN  TOE  TEST  DATA  FILE)" 

180  PRINT  TAB(S );"*••**•••••••• . .................. . . . 


190  INPUT  "<CR>  TO  CCNTINE"  ;XJ 

200  PRINT  TAB(J);"* . . . . . 

210  PRINT  TAB ( 5 ) ;  "TOE  PROGRAM  TAKES  TOE  DATA  FILE  CREATED  FROM  TOE  PDPll/04" 

220  PRINT  TAB( 5 ) ; "DISC.  AS  FILE  bmw. OUT;#,  AM)  READS  EACH  LINE  TO  DETERMINE" 
230  PRINT  TAB ( 5 ) ; "TOE  Z9  VALUE.  WEN  Z9-3  DATA  IS  EXTRACTED  FROM  THE  FILE  AND" 
240  PRINT  TAB ( 5 ) ;  "UNLOADING  COMPLIANCE  CALIBRATION  IS  USED  TO  MEASURE  DELTA  A' 
250  PRINT  TAB ( 5 ) ;  "THE  J  VALUES  ARE  CALCULATED  FROM  TOE  DATA  STORED  IN  TOE  FILE" 

260  PRINT  TAB ( 5  ) ; . . . . . . . . . 

270  INPUT  "<CR>  TO  CONTINUE*  ;X5  r 

280  PRINT  TAB(5 ) ; . . . . . . 

290  PRINT  TAB(5) ; 'SPECIMENS  FOR  VHICH  UNLOADING  COMPLIANCE  CALS.  ARE  AVAILABLE" 
300  PRINT  TAB(5) ; "SPECIMEN  GECMTRY  AND  CALIBRATION" ;TAB(55); "TYPE" 

310  PRINT  TAB(S)  ;*  MIL  J03  /  PAFBC  1984";TAB(S7);"1"  /  PRINT 
320  PRINT  TAB ( 5 ) ; '  IfltL  J03  /  EXPERftRNT"  ;TAB( 57 ) ;  "2"  /  PRINT 


290  PRINT  TAB(5) ; 'SPECIMENS  FOR  WICH  UNLOADING  COM 
300  PRINT  TAB(5) ; 'SPECIMEN  GECMTRY  AND  CALIBRATION" 
310  PRINT  TAB(S)  ;*  MIL  J03  /  PAFBC  1984*;TAB(S7);"1 
320  PRINT  TAB ( 5 ) ; '  IfltL  J03  /  EXPERIMENT" ;TAB(S7);"2 

380  PRINT  TAB(5)  ;"••*•••••• . . . •*•••••• 

390  INPUT  "INPUT  TOE  SPECIMEN  TYPE  ";TYPE%  /  PRINT 
400  IF  TYPESK2  GOTO  430 

410  PRINT  "INVALID  SPBCIVRN  TYPE  *  /  PRINT 
420  GOTO  280 


430  ON  TYPER  GOSUB  1860.1910 

440  PRINT  "POLYNOMIAL  OOEFFS  FOR  TYPE  "; TYPER;"  SPECIMEN" 
450  PRINT  NOTES 

460  PRINT  "X0-  '  ;X0  /  PRINT  "XI-  ";X1  /  PRINT  "X2-  ";X2 

470  PRINT  "X3-  ";X3  /  PRINT  "X4-  ";X4  /  PRINT 

480  INPUT  "IS  THIS  CALIBRATION  TOE  ONE  REQUIRED?  Y/N  ";A$ 

490  IF  AS-"N"  GOTO  280 

500  PRINT  /  PRINT 


510  PRINT  ••••••••• . .  STARTING  TO  RERUN  TOE  DATA 

520  PRINT  /  PRINT 

530  MAP  (A)  STRING  NN  -  132R 


540  INPUT  "WAT  IS  TOE  OUTPUT  DATA  FILE  NAIC,  VIZ  name  ";TS 

550  OPEN  TJ‘'.OUT"  FOR  INPUT  AS  FILE  *1R.  SEQUENT  I  AL.RECORDTYPE  ANY  A.  RECORDS  I ZE  132%,  MAP  A 
560  INPUT  "HW  MANY  READINGS  IN  TOE  FILE"  ;QR 

570  INPUT  "WAT  IS  THE  MERKLE-CORTEN  OOEFF  FOR  THIS  SPECIMEN  ";MX 
580  DIM  P(50%)  ,L(50%)  ,A(50%)  ,N( 50% ) 

590  DIM  R5(50%),DA5(50%),A5(50%).U5(50%).J5(50%) 

600  DIM  F5(50»).VS(50%).C5(S0%) 

610  Z9%-2%  /  N7%-0% 


620  GET  *1% 

630  SNJ-SEGJ(FN,43.5l) 
640  PRINT  SNJ ; "  BITS/KN" 
650  RNS-SEGS(N4,57,6l) 
660  PRINT  RNS ; "  BITS/lM" 


I 


670  TWJ-SEG5(W,  1 .13) 

680  TNS “TRMS  (TTWS ) 

690  PRINT  "SPECIVEN  THICKNESS  AND  WIDIH  SHOULD  BE  MEASURED  WITH  A  MICROMETER.  ” 
700  INPUT  "SPECIMEN  THICKNESS  (METRES)" ;B 
710  INS-SEGS(N4,14.19) 

720  INPUT  "SPECIMEN  WIDTH  (METRES  )";W 
.  730  VNS-SEGS (N4, 28 , 37 ) 

740  PRINT  "YOUNG’S  MXULAS-';VNJ;’P*’ 

750  M-VAL(VNS) 

760  INPUT  "DO  YOU  WISH  TO  CHANCE  YONG’S  MODULUS,  Y/N  ";SS 
770  IF.SS-"N"  GO  TO  800 

780  INPUT  "WHAT  IS  THE  NEW  YOUNG’S  MODULUS ,  . . .  IN  GP*  "  ;M 
790  M-M"1E*09 
800  PRINT  B.W.M 
810  NJ-VAL(SNJ) 

820  N4%“VAL(RNS ) 

830  K-N3/(N4%*1.0E*06) 

840  FOR  I*-l%  TO  Q % 

850  GET  #1% 

860  SNS-SEGS  (fN,  58,58) 

870  a-VAL(SNS)  r  ' 

880  IF  Z*-l%  GO  TO  920  ( 

890  IF  »-2%  THEN  N%-0% 

900  IF  ZX-(Z9%-1%)  THEN  GOSUB  1200 
910  IF  ZX-3%  THEN  GOSUB  950 
920  NEXT  I* 

930  IF  »-3%  THEN  GOTO  1200 
940  GO  TO  1620 

950  REM  DATA  ACQUISITION  FROM  PDP11  DATA  FILE  VF1 
960  IF  NKO0%  THEN  GOTO  1000 
970  PRINT  PRINT  PRINT 

980  PRINT  "UNLOADING  COMPLIANCE  DATA"  PRINT 
990  PRINT  "READING", "LOAD", "DISPL.", "AREA" 

1000  NlS"SEGS(fN.  1,5 ) 

1010  N2S-TRMS(N1S) 

1020  N6%-VAL(N2S) 

1030  P1S-SEGS (fW.  10,23) 

1040  P2J-TRMS(P1J) 

1050  P-VAL(P2S) 

1060  L1S-SBGS(M4.23.35) 

1070  L2S-TRMS(L1J) 

1080  L-VAL(L2S) 

1090  A1  J”SEGJ  (fN,  38,53) 

1100  A2S-TRMS  (AlS) 

1110  A-VAL(A2S) 

1120  P(NR)-P 
1130  L(N»)-L 
1140  A(N%)-A 
1150  N(N»)-N6% 

1160  PRINT  N(N*),P/N3.L/N4%,A/(N3*N4%) 

1170  N»-N»*l* 

1180  Z9%-»  /  N5%-N% 

1190  RETURN 

1200  PRINT  /  PRINT  /  PRINT 

1210  PRINT  "UNLOADING  COMPLIANCE  RERUN" 

1220  RBI  SETTING  SUMMATIONS  TO  ZERO 
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1230  El-0  /  E2-0  /  E3-0  /  E4-0  /  E5-0  /  E6-0 

1240  E7-0  /  E8-0  /  E9-0  /  E10-0  /  F%-0 

1230  PRINT  'MAX.  NO.  OF  UNLOADING  POINTS  IS  ";N5% 

1260  INPUT  "HOW  MANY  OFF  THE  MAX.  LOAD  ARE  TO  BE  DISCARDED  ";0* 

1270  IF  OX>(NS%-10%)  THEN  PRINT  'ERROR  IN  INPUT,  TRY  AGAIN"  /  GOTO  1250 

12*0  IF  O%-0  THEN  PRINT  "DISCARD  OF  ZERO  POINTS  IS  OK?  <Y/N>  *  /  INPUT  ZS  /  IF  ZSO"Y"  GOTO 
1290  FOR  C*-(0VO*)  TO  (N5%-1»)  STEP  M% 
i:-.  1300  El-  El *P(C%)  /  E2-E2*L(C»)  /  E3-E3*P(CX)*L(C*) 

W  1310  E4-E4*L(CX)  2»  /  E6-E6*P(CX)  2% 

H  1320  F%-F%*  1* 

:  1330  NEXT  CX 

1340  REM  CALCULATION  OF  SLOPE 

1  1350  ES-E3-((E1*E2)/F»)  /  E9-E4-((E2  2%)/F%) 

1360  E5-E8/E9 

1370  REM  CALCULATION  OF  CORRELATION  OOEF. 

1380  E10-E6-((E1  2*)/F%) 


1390  E7-(E8  2)/(E9*E10) 

1400  V- (K*M*B )/ES  /  REM  NORMALIZED  COMPLIANCE 

1410  V5(N7%)-V  /  C5(N7%)-E7 

1420  R-XO ♦XI  *V*X2 *V  2»*X3*V  3%*X4*V  4% 

1430  IF  N7%>-1%  GOTO  1450  ^ 

1440  Y1-R*W  /  Y0-Y1  /  GOTO  1460  1 

1450  Y1-R*W 

1460  PRINT  "CHECK  OF  NO  OF  UNLOADING  POINTS  ";F%;*  SHOULD  BE  ";N5X-OK 
1470  PRINT  'CRACK  LENGTH-  ";Y1 ;"  METRES" 

1480  PRINT  "READING", "DATA  POINTS", "NORM.  CCMPL.  " , "CDRR.OOEF.  " ."AREA" , ”A/W , "DELTA  A","J" 

1 490  J-MOC*A(CK )/(B* (W-Yl ) *N3 *N4%*  1 000 ) 

1500  UX-N6X-NSV1X 

1510  PRINT  L».F».Y.E7.A(CJI)/(N3*N4»),R.Y1-Y0.J 
1520  Z9%-» 

1530  U5(N7%)-U*  /  A5(N7X)-A(ai)/(N3*N4%)  /  R5(N7»)-R 
1540  DA5(N7»)-Y1-Y0  /  J5(N7%)-J  /  F5(N7»)-F% 

1550  INPUT  *<CR>  TO  PROCEED  TO  NEXT  UNLOADING  SET  ";XS 
1560  N7%-N7%*  1% 

1570  IF  IVXJ*  THEN  GOTO  1610 

1580  PRINT  "LAST  UNLOADING  DATA  SET" 

1590  INPUT  "SET  PRINTER  AM)  <CR>  FOR  PRINT  OUT  OF  UNLOADING  COMPLIANCE  CALCS'; XI 
1600  GOTO  1620 
1610  RETURN 
1620  CLOSE  *1% 

1630  PRINT  /  PRINT  /  PRINT 

1640  PRINT  "SUMMARY  OF  UNLOADING  COMPLIANCE  CALCULATIONS" 

1650  PRINT  "RUN  USING  JRERLN.BAS  ON  THE  VAX1 1/780" 

1660  PRINT  "USING  PDPll/04  DATA  FILE  ";TS;"(RUN  THROUGH  V.BAS)" 

1670  PRINT 

1680  PRINT  "YOUNG'S  MODULUS  FOR  THESE  CALCULATIONS  -";M;"Pa" 

1690  PRINT  "SPECIMEN  THIONESS-";B;"metrej" 

1700  PRINT  "SPECIMEN  WIDTH  -"  ;W;  "me  t  re  j"  /  PRINT 
1710  PRINT  NOTES 

1720  PRINT  "POLYNOMIAL  FIT  COEFFS" 

1730  PRINT  "X0-  *;X0;"  XI-  ";X1 ;"  X2-  ";X2 
1740  PRINT  "X3-  ";X3;"  X4-  ";X4 
1750  PRINT 

1760  PRINT  "MERKLE-CORTEN  COEFF-  ";MX 
1770  PRINT  /  PRINT 

1780  PRINT  "READING",  "DATA  POINTS' ."NORM  CCMPL",  "OORR  OOEF .  "  ,  "AREA"  .  "A/W"  ,  "DELTA  A"  , "  J  VALUE 
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1790 

- 

1800 

1810 

• 

1820 

1830 

1840 

M  • 

18S0 

1860 

1870 

1880 

1890 

1900 

l 

1910 

1920 

1930 

1940 

1950 

*  #-.»i 

*•  • 

1960 

2110 

*: ,  - 

2120 

2130 

•- 

2140 

2150 

2160 

2170 

2180 

V  4 

2190 

&  * 

2200 

2210 

2220 

2230 

2240 

2250 

i 

2260 

yi 

2270 

2280 

2290 

2300 

•  ’  -4 

2310 

2320 

2330 

2340 

2350 

FOR  NX-OX  TO  N7X- 1% 

PRINT  US  (NX ) .  F5  (NX ) , V5  (NX )  ,C5  (NX) . A5  (NX ) ,  R5  (NX )  ,QA5  (NX ) .  J  J  (NX ) 

NEXT  NX 
PRINT  /  PRINT 

PRINT  THIS  FILE  WILL  BE  STORED  AS  *;TS;".DAT" 

PRINT  TON  J  PRINT  TO  RETRIEVE  THIS  FILE" 

GOTO  1110 

REM  COMPLIANCE  CALIBRATIONS  (POLYNOMIAL  FITS) 

NOTES -'JO 3  CAL.  FROM  PAFBC  19*4" 

X0-7.4857727E-02  /  Xl-1 . 7894292E-02  /  X2— 2. 3063905E-04 

X3-1.5006292E-06  /  X4--3. 7869343E-09 

RETURN 

REM 

NOTES -"JO 3  CAL.  FROM  CLARK  EXPT.  DATA" 

XO-1 . 0697574E-01  /  Xl-l . 47S4694E-02  /  X2--1 . 6242350E-04 

X3-8 . 7544527E-07  /  X4— 1 . 7795666E-09 

RETURN 

REM 

REM  PROQtAM  TO  WRITE  A  SUMMARY  OF  THE  JRERLN  CALCULATIONS  AS  A  FILE  name.DAT 
OPEN  TS-" .DAT"  FOR  OUTPUT  AS  FILE  *2%,  ORGANIZATION  SEQUENTIAL,  RECORETIYPE  ANY 
PRINT  #2%,  TS 

PRINT  *2X,M  *' 

PRINT  #2%,B 
PRINT  #2%,W 
PRINT  #2*  .NOTES 

PRINT  #2%.X0  /  PRINT  *2»,X1  /  PRINT  #2%,X2  /  PRINT  #2X.X3  /  PRINT  #2».X4 
PRINT  #2*  .MX 
RX-N7X- 1% 

PRINT  #2*,R% 

FOR  NX-0%  TO  R% 

PRINT  #2%. US  (NX) 

PRINT  #2%.F5(NX) 

PRINT  #2X,Y5(NX) 

PRINT  *2%.CS(NX) 

PRINT  #2X,AJ(NX) 

PRINT  *2X,R5(NX) 

PRINT  #2% , DAS  (NX ) 

PRINT  #2% ,  J 5  (NX ) 

NEXT  NX 
CLOSE  #2X 
PRINT  /  PRINT 

PRINT  "END  OF  JRERLN. BAS . program  by  D.S.  Saunders  and  I.A.Bnrch  1985" 

END 
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(3)  ACCCSS1 . BAS 


REM  .  .  .AUSTRALIAN  GOVERFMENT  DEFENCE  SCIENTIFIC  AND  TECHNOLOGY-  ORGINI2ATICN 
110  REM  ...MATERIALS  RESEARCH  LABORATORIES 

120  REM  THIS  PROGRAM  WILL  ACCESS  OUTPUT  FILES  FOR  TEST  DATA  FOR  THE  VAX  1 1/780 

130  PRINT  TAB (3  );"••••*•••••••••• . . . . . . . 

140  PRINT  TAB(5) -.'ACCESS  1.  A  PROOtAM  TO  CONVERT  A  J  INTEGRAL  TEST  DATA  FILE' 

150  PRINT  TAB ( 5 ) ; ' INTO  A  PLOTTING  DATA  FILE  FOR  USE  WITH  THE  KENNETT  PLOT  ROUTINE" 

160  PRINT  TAB(S) . . . * . * . 

170  PRINT  TAB ( 5 ) ;  "THE  PROGRAM  TAKES  THE  DATA  FILE  CREATED  FROM  THE  PDPll/04" 

180  PRINT  TAB(S ) ; "DISK,  AS  FILE  SSSS.OUT;#,  AND  WRITES  A  NEW  FILE  ( FORCE/DI S PL.  )" 

190  PRINT  TAB ( 5 ) ; "NAME  THIS  FILE  SSSS.GRA;#" 

200  PRINT  TAB ( 5 ) ; "THE  PROGRAM  WILL  RECORD  THE  VAX  DATA  FILE  FROM  WHICH  SSSS.GRA;#  IS  CREATED” 

210  PRINT  TAB(5 ) ;  . . . . . . •*" 

220  MAP  (A)  STRING  NN  -  132% 

230  INPUT  "WHAT  IS  THE  OUTPUT  FILE  NAME.  VIZ  SSSS.OUT;#.  ";TJ 

240  OPEN  T$  FOR  INPUT  AS  FILE  #1%,  SEQUENTIAL, RECORDTVPE  ANY  &, RECORDS  I ZE  132%,  MAP  A 
250  INPUT  "WHAT  IS  THE  PLOTTING  FILE  NAME.  VIZ  SSSS.GRA;#.  ";SS 
260  OPEN  SS  FOR  OUTPUT  AS  FILE  #2%, ORGANIZATION  SEQUENTIAL  VARIABLE 
270  INPUT  "WHAT  IS  THE  FIRST  READING  NUMBER  YOl/  WANT  TO  ACCESS"  ;P% 

280  INPUT  "WHAT  IS  THE  LAST  READING  NUMBER  YOU  WANT  TO  ACCESS'  ;Q% 

290  GET  *1% 

300  SNS-SEGS(M4,43 ,51 ) 

310  PRINT  SNS 

320  RNS-SEGS(NN, 57,61) 

330  PRINT  RNS 
340  N3 “VAL ( SNS ) 

350  N4%-VAL%(RNS) 

360  PRINT  TAB  (20);  "SELECTED  FORCE/D  IS  PLACEMENT  DATA  FRCM";TS 
370  PRINT  TAB(20); "WRITTEN  INTO  PLOTTING  FILE  ";SS 
380  PRINT  TAB (20); "SCALING  CONSTANTS  FOR  THE  TEST" 

390  PRINT  TAB(20);"FOR  LOAD,  RO-  "  ;N3  ;TAB(50 ) ;  "BITS/KN" 

400  PRINT  TAB (20); "FOR  DISPL,  Rl-  " ;N4%;TAB (50 ) ; "BITS/VM" 

410  PRINT  /  PRINT 

420  PRINT  #2%. "DISPLACEMENT  IN  VM" 

430  PRINT  #2% ,  1% ; " .  " 

440  PRINT  #2%. "FORCE  IN  KN" 

450  PRINT  #2% , 1% ; " . " 

460  US-LEFTS (TS,4) 

470  PRINT  #2%. "DATA  FRCM  ORIGINAL  CN  PDP1 1/04  DISK.  AS  ";US;".DAT" 

480  PRINT  TAB(5); "READINGS  NO  ";P%;"  TO  ";Q%;  "INCLUSIVE" 

490  PRINT  /  PRINT 
500  B%—  ((Q*-P%)M) 

510  PRINT  *2%,B% 

520  PRINT  #2%,  "PLOTTING  DATA  FROM  ";TS ;  "READING  NOS  ";P%;"TO  ";Q%;"INCL" 

530  FOR  1%  -  1%  TO  2000% 

540  GET  #1% 

550  NNWS-RKHTS(fW,7) 

560  N1S-LEFTS(NNWS,17) 

570  N2S-SEGS(NWS.18,33) 

580  Nl-VAL(NIS) 

590  N2-VAL(N2S) 

600  IF  I%<P%  THEN  640 
610  IF  I%>Q*  THEN  650 
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620  PRINT  N2/N4%,N1/N3 

630  PRINT  #2%,N2/N4*;”,’;N1/N3 

640  NEXT  1% 

630  CLOSE  #1%  /  CLOSE  #2% 

660  PRINT  ’EM)  OF  ACCESS. BAS. . 
670  DO 


Program  by  D.S.Saunder 


APPENDIX  2 


Ustlngs  of  Elastic  Unloading  Compliance  for  a  Range 
of  Specimen  Geometries 


MKL  J03  (compact  Tension  Specimen) 
(1)  PAFEC  -  generated  (1984)  Ill 


VLL 

EB  — 

a/W 

18.0958 

0.3 

19.0547 

0.32 

2 0.3038 

0.34 

21.7481 

0^36 

23.563  i 

0.38 

25.4099 

0.4 

27.4611 

0.42 

29.7407 

0.44 

32.4349 

0.46 

35.3034 

0.48 

38.5184 

0.5 

42.1361 

0.52 

46.2305 

0.54 

50.2362 

0.58 

62.4079 

0.6 

69.5932 

0.62 

78.0284 

0.64 

88.0301 

0.66 

100.013 

0.68 

113.907 

0.7 

131.604 

0.72 

153.631 

0.74 

181.487 

0.76 

217.361 

0.78 

264.575 

0.8 

(2)  PAFEC-generated  (1985)  (2),  Point  Pin  Loading 


0.52 

0.56 

0.60 

0.64 

0.68 

0.72 

0.76 


The  second  PAFEC-generated  compliance  calibration  [2]  is  now 
considered  to  be  more  accurated  than  (l),  above.  For  this  calibration  the 
original  mesh  was  modified  slightly  and  the  normalised  compliance  data  is 
close  to  the  experimental  data  given  in  [31  and  reproduced  below. 


44.74 

54.43 

67.30 

85.00 

109.9 

147.5 

207.3 


__  LL 

EB  - 

P 

a/w  , 

24.78 

0.378 

26.64 

0v^393 

29.58  *’ 

0.419 

40.76 

0.493 

41.6 

0.495 

41.65 

0.497 

40.58 

0.498 

41.63 

0.498 

42.76 

0.507 

44.5 

0.517 

53.47 

0.544 

55.15 

0.561 

57.75 

0.565 

59.43 

0.567 

90.6 

0.644 

110.3 

0.670 

154.0 

0.722 

167.0 

0.735 

209.0 

0.763 

212.1 

0.764 

215.9 

0.767 

225.7 

0.768 

224.9 

0.773 

242.5 

0 . 778 

251.9 

0.784 

278.0 

0.793 

(B) 


H  Tension  Specimen) , 
V 

141,15) 

„  LL 

EB  — 

a/w 

8.6 

0.2 

11.18 

0.25 

14.28 

0.30 

18.09 

0.35 

A2.2 


22.86 

0.40 

28.96 

0.45 

36.99 

0.50 

47.90 

0.55 

63.35 

O 

VO 

• 

0 

86.36 

0.65 

122.80 

0.70 

185.40 

0.75 

304.6 

0.80 

578.36 

0.85 

1389.07 

0.90 

5952.20 

0.95 
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Lowing 


A2.3 


A1 


Area 


R 


S  * 


2 


A3 


E5 


V 

R 


Y1 

D5 

P5 

D9 


This  is  the  area  under  the 
linear-elastic  part  of  the 
load/ load-line-displacement 
record  and  is  calculated  as  the 
area  of  the  P/D  triangle  to  P7. 


Area  This  area  is  continuously  up¬ 

dated  for  each  displacement 
increment,  viz  Z9-1 ,  Z9-2,  but 
not  when  Z9-3  or  Z9-4.  A3  at 
the  final  load  increment  before 
the  unload  is  written  into  the 
virtual  files  VFl,  VF2.  A3  is 
not  increased  again  until  Z9 
returns  to  a  value  of  2. 


Raw  compliance  This  value  from  the  linear 

regression  of  the  unloading 
data. 


Normalized  compliance  ^E5  xExBxPxr— 

RO 

a/w  This  value  is  calculated  from 

the  unloading  compliance 
calibration  polynomial  (f(v)). 


Crack  length 
Displacement  (bits) 
Load  (bits) 
Displacement  (bits) 


R  *  w 

Average  of  10  values 

Average  of  10  values 

Displacement  at  the  commencement 
of  unloading. 


(l)  See  listing  3 
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Example  of  an  Output  from  JRERUN.BAS 


